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Abstract 

Consider a sub-Riemannian geometry (U, D,g) where U is a neighbor- 
hood at in 1R™, D is a rank-2 smooth (C°° or C w ) distribution and g is 
a smooth metric on D. The objective of this article is to explain the role 
of abnormal minimizers in SR-geometry. It is based on the analysis of the 
Martinet SR-geometry. 
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Resume On considere un probleme sous-Riemannien (U,D,g) ou U est un 
voisinage de dans IR" , D une distribution lisse de rang 2 et g une metrique lisse 
sur D. L'objectif de cet article est d'expliquer le role des geodesiques anormales 
minimisantes en geometrie SR. Cette analyse est fondee sur le modele SR de 
Martinet. 
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1 Introduction 

Consider a smooth control system on H™ : 

q(t) - f(q(t),u(t)) (1) 

where the set of admissible controls hi is an open set of bounded measurable 
mappings u defined on [0, T(u)] and taking their values in H m . We fix q(0) = q 
and T(u) = T and we consider the end-point mapping E : u £ hi i — > q(T, qo, u), 
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where q(t,qo,u) is the solution of (JIJ associated tougW and starting from go 
at t = 0. Wc endow the set of controls defined on [0, T] with the L°°-topology. 
A trajectory q(t, q$, u) denoted in short q is said to be singular or abnormal on 
[0, T] if u is a singular point of the end-point mapping, i.e, the Frechet derivative 
of E is not surjective at u. 

Consider now the optimal control problem : min J„ T f°{q(t), u{t))dt, where 

f is smooth and q(t) is a trajectory of JQ) subject to boundary conditions : 
q(0) G Mo and q(T) £ Mi, where Mq and M\ are smooth submanifolds of 
R™. According to the weak maximum principle |37| . minimizing trajectories are 
among the singular trajectories of the end-point mapping of the extended system 
in H n+1 : 

q(t) = f(q(t),u(t)) (2) 
q°(t) = f(q(t),u(t)) 

and they are solutions of the following equations : 
. dH v . dH v dH v 

q = ^—, p = — ^— , -5— = o 3 

op oq ou 

where H v =< p, f(q,u) > +vf°(q,u) is the pseudo-Hamiltonian, p is the ad- 
joint vector, <, > the standard inner product in H™ and v is a constant which 
can be normalized to or —1/2. Abnormal trajectories correspond to v = ; 
their role in the optimal control problem has to be analyzed. Their geometric 
interpretation is clear : if C denotes the set of curves solutions of Q, they 
correspond to singularities of this set and the analysis of those singularities is 
a preliminary step in any minimization problem. This problem was already 
known in the classical calculus of variations, see for instance the discussion in 
|H] and was a major problem for post-second war development of this disci- 
pline whose modern name is optimal control. The main result, concerning the 
analysis of those singularities in a generic context and for affine systems where 

m 

f(q,u) = F (q) + ^\nFi(q) are given in ^21 an d H3- The consequence of this 

t=i 

analysis is to get rigidity results about singular trajectories when m = 1, that is 
under generic conditions a singular trajectory 7 joining qo = 7(0) to q\ = j(T) 
is the only trajectory contained in a C°-neighborhood of 7 joining q to qi in 
time T (and thus is minimizing). 

In optimal control the main concept is the value function S defined as follows. 
If go j <7i and T are fixed and 7 is a minimizer associated to u 7 and joining q$ to 
qi in time T, we set : 

S(q ,q u T)= [ / ( 7 (i),« 7 (i))* 



The value function is solution of Hamilton- Jacobi-Bellman equation and one of 
the main questions in optimal control is to understand the role of abnormal 
trajectories on the singularities of S. 



2 



The objective of this article is to make this analysis in local sub-Riemannian 
geometry associated to the following optimal control problem : 

pT m 

min / Ui(t)dt , u = (u u . . . , u m ) 

u{.) Jo ^ 

subject to the constraints : 

tn 

g(t)=£« i (t)F i (g(*)) (4) 
»=i 

q G U C R n , where {Fi,--- , F m } are m linearly independant vector fields 
generating a distribution D and the metric g is defined on D by taking the F{ 's 
as orthonormal vector fields. The length of a curve q solution of Q on [0, T] 

and associated to u e W is given by : L(q) = I \J2 u i(t)) dt and the SR- 

distance between go and q\ is the minimum of the length of the curves q joining 
qo to q\. The sphere >S'(goj r ) with radius r is the set of points at SR-distance r 
from go- If & nv pairs qo, q% can be joined by a minimizcr, the sphere is made of 
end-points of minimizers with length r. It is a level set of the value function. 

It is well known (see pQ) that in SR-geometry the sphere S(qo, r) with small 
radius r has singularities. For instance they are described in 0] in the generic 
contact situation in M 3 and they are semi-analytic. Our aim is to give a geomet- 
ric framework to analyze the singularities of the sphere in the abnormal direc- 
tions and to compute asymptotics of the distance in those directions. We analyze 
mainly the Martinet case extending preliminary calculations from |21 EH ■ The 
calculations are intricate because the singularities are not in the subanalytic 
category even if the distribution and the metrics are analytic. Moreover they 
are related to similar computations to evaluate Poincare return mappings in the 
Hilbert's 16th problem (see |35l I38p using singular perturbation techniques. 

The organization and the contribution of this article is the following. 

In Sections 2 and 3 we introduce the required concepts and recall some 
known results from [51[7|E] to make this article self-contained. 

In Section 2, we compute the singular trajectories for single-input affinc 
control systems : q = Fo(q)+uFi(q), using the Hamiltonian formalism. Then we 
evaluate under generic conditions the accessibility set near a singular trajectory 
to get rigidity results and to clarify their optimality status in SR geometry. 

In Section 3 we present some generalities concerning SR geometry. 

In Section 4 we analyze the role of abnormal geodesies in SR Martinet geom- 
etry. We study the behaviour of the geodesies starting from in the abnormal 
direction by taking their successive intersections with the Martinet surface filled 
by abnormal trajectories. This defines a return mapping. To make precise com- 
putations we use a gradated form of order where the Martinet distribution is 

2 

identified to Ker u), u> = dz — ^-dx, the Martinet surface to y — 0, the abnor- 
mal geodesic starting from to t i— > (t, 0,0) and the metric is truncated at : 
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g = (1 + ay) 2 dx 2 + (1 + fix + ■yy) 2 dy 2 , where a, (3, 7 are real parameters. The 
geodesies equations project onto the planar foliation : 



6" + sm6 + e/3cos8 9' + e 2 asm8(acos9 - (3sm6) = 



(5) 



where e — -A= is a small parameter near the abnormal direction which projects 
onto the singular points 6 = kir. Here the Martinet plane y = is projected 
onto a section E given by : 



Equation (JjJ represents a perturbed pendulum ; to evaluate the trace of 
the sphere S(0,r) with the Martinet plane we compute the return mapping 
associated to the section E. Near the abnormal direction the computations are 
localized to the geodesies projecting near the separatrices of the pendulum. 

In order to estimate the asymptotics of the sphere in the abnormal direction 
we use the techniques developped to compute the asymptotic expansion of the 
Poincare return mapping for a one-parameter family of planar vector fields. This 
allows to estimate the number of limit cycles in the Hilbert's 16th problem, see 
|38) . Our computations split into two parts : 

1. A computation where we estimate the return mapping near a saddle point 
of the pendulum and which corresponds to geodesies close to the abnormal 
minimizer in C 1 -topology. 

2. A global computation where we estimate the return mapping along geodesies 
visiting the two saddle points and which corresponds to geodesies close to 
the abnormal minimizer in C° -topology, but not in C 1 -topology (see |42| 
for a general statement). 

Our results are the following. 

• If P = 0, the pendulum is integrable and we prove that the sphere belongs 
to the log-exp category introduced in [T§|, and we compute the asymptotics 
of the sphere in the abnormal direction. 

• If (3 7^ 0, we compute the asymptotics corresponding to geodesies C 1 close 
to the abnormal one. 

We end this Section by conjecturing the cut-locus in the generic Martinet 
sphere using the Liu-Sussmann example |29| . 

The aim of Section 5 is to extend our previous results to the general case 
and to describe a Martinet sector in the n-dimensional SR sphere. 

First of all in the Martinet case the exponential mapping is not proper and 
the sphere is tangent to the abnormal direction. We prove that this property is 
still valid if the abnormal minimizer is strict and if the sphere is C^-stratifiable. 

Then we complete the analysis of SR geometry corresponding to stable 2- 
dimensional distributions in H 3 by analyzing the so-called tangential case. We 



6' = e (a cos9 + (3 am 9) 



(6) 
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compute the geodesies and make some numerical simulations and remarks about 
the SR spheres. 

The flat Martinet case can be lifted into the Engel case which is a left- 
invariant problem of a 4-dimensional Lie group. We give an uniform parametriza- 
tion of the geodesies using the Weierstrass function. Both Heisenberg case and 
Martinet flat case can be imbedded in the Engel case. 

The main contribution of Section 5 is to describe a Martinet sector in the 
SR sphere in any dimension using the computations of Section 4. We use the 
Hamiltonian formalism (Lagrangian manifolds) and microlocal analysis. This 
leads to a stratification of the Hamilton- J acobi equation viewed in the cotangent 
bundle. 
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2 Singular or abnormal trajectories 

2.1 Basic facts 

Consider the smooth control system : 

q(t) = f(q(t),u(t)) (7) 

and let q{t) be a trajectory defined on [0, T] and associated to a control u G U. 
If we set : 

A{t) = % (?(*)'«(*)) - B ® = |£ (?(*)»«(*)) 
the linear system : 

Sq(t) = A(t)Sq(t) + B(t)5u(t) (8) 

is called the linearized or variational system along (q,u). It is well known, see 
[T2] that the Frechet derivative in L°°-topology of the end-point mapping E is 
given by : 



E'{v) = $(T) [ ^- 1 (s)B(s)v(s)ds 
Jo 



where $ is the matricial solution of : $ = ^4$, with $(0) = id. 

Hence (q, u) is singular on [0, T] if and only if there exists a non-zero vector 
p orthogonal to Im E'(u), that is the linear system (jSJ) is not controllable on 

[o,n 

If we introduce the row vector : p(t) — p<l>(T)$ a standard compu- 

tation shows that the triple (q,p,u) is solution for almost all t £ [0,T] of the 
equations : 

. , . df df 

q = f(x,u), p = -p— , Pj- = 
oq ou 
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which takes the Hamiltonian form 



. _ aff . _ _dH_ dH_ _ (| 

^ dp ' P dq 1 du 

where H{q 1 p,u) =< p,f(q,u) > . The function iJ is called the Hamiltonian 
and p is called the adjoint vector. 

This is the parametrization of the singular trajectories using the maximum 
principle. We observe that for each t &]0,T[, the restriction of (q,u) is singular 
on [0, t] and at t the adjoint vector p(t) is orthogonal to the vector space K(t) 
image of L°°[Q, t] by the Frechet derivative of the end-point mapping evaluated 
on the restriction of u to [0, t]. The vector space K(t) corresponds to the first 
order Pontryagin's cone introduced in the proof of the maximum principle. If 
t e]0,T] we shall denote by k(t) the codimension of K(t) or in other words 
the codimension of the singularity. Using the terminology of the calculus of 
variations k(t) is called the order of abnormality. 

The parametrization by the maximum principle allows the computation of 
the singular trajectories. In this article we are concerned by systems of the 
form : 

q{t) = F {q{t))+u{t)F x (q(t)) (10) 
and the algorithm is the following. 

2.2 Determination of the singular trajectories 
2.2.1 The single input affine case 

It is convenient to use Hamiltonian formalism. Given any smooth function 
H on T*U, H will denote the Hamiltonian vector field defined by H. If Hi, 
Hi are two smooth functions, {H X ,H2} will denote their Poisson bracket : 
{Hi,H2\ = dH\(H2). If X is a smooth vector field on U, we set H =< p, X(q) > 
and H is the Hamiltonian lift of X. If X\,X2 are two vector fields with 
Hi =<p,Xi(q) >, i = 1,2 we have : {Hi, Hz} =< p, [Xi, X 2 ](q) > where the 

dX dX 

Lie bracket is : [Xi,X 2 ](g) = — — (q)X2(q) = — (q)X 1 (q). We shall denote 

oq oq 

by H =< p, F (q) > and H 1 =< p, Fi(q) > . 

If f(q, u) = Fo(q) + uFi(q), the equation © can be rewritten : 

8H dHi (dH Q dHi 

U— , P= — — r + U— 



op op V oq 

Hi=0 for all t € [0, T] 

We denote by z = (q,p) £ T*U and let (z,u) be a solution of the above 
equations. Using the chain rule and the constraint : Hi = 0, we get : 

0=4- Hi{z(t)) = dHi(z(t))H (z(t)) + u(t)dHi(z(t)) Hi(z(t)) for a.e. t 
dt 



6 



This implies : = {-Hi, Ho}(z(t)) for all t. Using the chain rule again we get : 

= {{H 1 ,H o },H o }(z(t)) + u(t){{H 1 ,H o },H 1 }(z(t)) for a.e. t 

This last relation enables us to compute u(t) in many cases and justifies the 
following definition : 

Definition 2.1. For any singular curve (z, u) : J — [0, T] i — > T*U xR, TZ(z, u) 
will denote the set {t e J, {{H Q , Hi}, Hi}(z(t)) ^ 0}. The set K(z,u) possibly 
empty is always an open subset of J. 

Definition 2.2. A singular trajectory (z, u) : J — > T*U x H is called of order 
two if TZ(z, u) is dense in J. 

The following Proposition is straightforward : 

Proposition 2.1. If{z,u) : J i — > T*MxlR is a singular trajectory andlZ(z,u) 
is not empty then : 

1. z restricted to 7Z(z,u) is smooth ; 

Conversely, let (Fo, Fi) be a pair a smooth vector fields such that the open 
subset fi of all z S T*{7 such that {{i/ , -ffiK^) 7^ is not empty. If 

\\H H 1 H 1 

iJ : £1 1 — > H is the function Ho + i?i then any trajectory of 

{{Hi,H },Hi\ 

H starting at t = from the set iJi = i/o} = is a singular trajectory of 
order 2. 

This algorithm allows us to compute the singular trajectories of minimal 
order. More generally we can extend this computation to the general case. 

Definition 2.3. For any multi-index a £ {0, 1}", a = [a\, . . . ,a n ) the func- 
tion H a is defined by induction by : H a = {H( ai H a „}. A singular 

trajectory (z,u) is said of order k ^ 2 if all the brackets of order m ^ k : Hp, 
with P = . . . , (3 m ), Pi = 1 are along z and there exists a — (1, a>2, ■ ■ ■ , Qffe) 
such that H a i(z) is not identically 0. 

The generic properties of singular trajectories are described by the following 
Theorems of |13) . 

Theorem 2.2. There exists an open dense subset G of pairs of vector fields 
(Fq,Fi) such that for any couple (Fq,Fi) £ G, the associated control system has 
only singular trajectories of minimal order 2. 

Theorem 2.3. There exists an open dense subset Gi in G such that for any 
couple (Fq, Fi) in G\, any singular trajectory has an order of abnormality equal 
to one, that is corresponds to a singularity of the end-point mapping of codi- 
mension one. 
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2.2.2 The case of rank two distributions 



Consider now a distribution D of rank 2. In SR-geometry we need to compute 
singular trajectories t i — ► q(t) of the distribution and it is not restrictive to 
assume the following : t i — ► is a smooth immersion. Then locally there 
exist two vector fields Fi,F 2 such than D = Span {Fi, F2} and moreover the 
trajectory can be reparametrized to satisfy the associated affine system : 

q(t) =u 1 (t)F 1 (q(t))+u 2 (t)F 2 (q(t)) 

where u\{t) = 1. It corresponds to the choice of a projective chart on the control 
domain. 

Now an important remark is the following. If we introduce the Hamiltonian 

2 

lifts : Hi =< p, Fi(q) > for i = 1, 2, and H — ^^UiHi the singular trajectories 

i=l 

are solutions of the equations : 

OH dH dH 

oq op ou 

dH 

Here the constraints — — = means : H\ = and H 2 = 0. This leads to the 
ou 

following definition : 

Definition 2.4. Consider the affine control system : q = F\ + uF 2 . A singular 
trajectory is said exceptional if it is contained on the level set : H = 0, where 
H =< p, Fi > +u < p,F 2 > is the Hamiltonian. 

Hence to compute the singular trajectories associated to a distribution we 
can apply locally the algorithm described in the affine case and keeping only 
the exceptional trajectories. An instant of reflexion shows that those of minimal 
order form a subset of codimension one in the set of all singular trajectories 
because H is constant along such a trajectory and the additional constraint 
Hi = has to be satisfied only at time t — 0. Hence we have : 

Proposition 2.4. The singular arcs of D are generically singular arcs of order 
2 of the associated affine system. They are exceptional and form a subset of 
codimension one in the set of all singular trajectories. 



2.3 Feedback equivalence 

Definition 2.5. Consider the class S of smooth control systems of the form : 

q(t) = f(q(t),u(t)), ?eE", ueR m 

Two systems f(x,u) and f'{y,v) are called feedback equivalent if there exists a 
smooth diffeomorphism of the form : $ : (x, u) 1 — ► (y, v), y = <p(x), v = ip(x, u) 
which transforms / into /' : 

d$(x)f{x,u) = f'(y,v) 
and we use the notation /' = $ * /. 
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Here we gave a global definition but there are local associated concepts which 
are : 

• local feedback equivalence at a point (xo, uq) <E H" x lR m . 

• local feedback equivalence at a point Xq of the state-space. 

• local feedback equivalence along a given trajectory : q(t) or (q(t),u(t)) of 
the system. 

This induces a group transformation structure called the feedback group G / 
on the set of such diffeomorphisms. For affine systems we consider a sub-group 
of Gf which stabilizes the class. This leads to the following definition. 

Definition 2.6. Consider the class of m-inputs affine control systems : 
J (t) = F (q(t) + F(q(t))u(t) 

m 

where F(q)u — J>UiFj(q). It is identified to the set A m +i — {F , Fi,..., F m } 

i-l 

of (m + l)-uplets of vector fields. The vector field F is called the drift. Let 
D be the distribution defined by D — Span {Fi(q), . . . , F m (q)}. We restrict the 
feedback transformations to diffeomorphisms of the form $ = (<p(q), ip(q,u) = 
a(q) + j3(q)u), preserving the class A. We denote by G the set of triples (f, a, (3) 
endowed with the group structure induced by Gf. 

We observe the following : take (F ,F) g A and $ = (ip,a,f3) £ G, then 
the image of (F , F) by <!> is the affine system (Fq, F') given by : 

(i) F^ cp * {F + R0) 

(ii) F' = Lp * F.(3. 

In particular the second action corresponds to the equivalence of the two 
distributions D and D' associated to the respective systems. 
The proof of the following result is straightforward, see [S] . 

Proposition 2.5. The singular trajectories are feedback invariants. 

Less trivial is the assertion that for generic systems, singular trajectories will 
allow to compute a complete set of invariants, see jS] for such a discussion. 

2.4 Local classification of rank 2 generic distributions D 
in JR 

We recall the generic classification of rank 2 distributions in M 3 , see 05] i with 
its interpretation using singular trajectories. Hence we consider a system : 

q(t) =u 1 (t)F 1 (q(t))+u 2 (t)F 2 (q(t)) 
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q = (x, y, z). We set D = Span {Fx,F 2 } and we assume that D is of rank 2. Our 
classification is localized near a point qo € R 3 and we can assume qo — 0. We 
deal only with generic situations, that is all the cases of codimension ^ 3. We 
have three situations which can be distinguished using the singular trajectories. 

Introducing Hi =< p, Fi(q) >,i = 1,2, a singular trajectory z = (q,p) must 
satisfy : 

H 1 = H 2 = {Hi,H 2 } = 

and hence they are contained in the set M : {q G R 3 ; det (F\, F 2 , [Fx, F 2 ]) = 0} 
called the Martinet surface. The singular controls of order 2 satisfy : 

Ul {{Hx,H 2 }, Hx} + u 2 {{H u H 2 }, H 2 } = 

We define the singular set S = SxDS 2 where 5, = {det(F x , F 2 , [[Fx, F 2 ], F z }) = 0}. 
We have the following situations. 

Case 1. Take a point qo £ M, then through q there passes no singular arc. 
In this case D is (C°° or C u )- isomorphic to Ker a, with a — ydx + dz. For this 

d 

normalization da = dy A dx (Darboux) and — is the characteristic direction. 
This case is called the contact case. 

Case 2. (Codimension one). We take a point q £ M\S. Since q $ S, we 
observe that M is near qo a smooth surface. This surface is foliated by the 
singular trajectories. A smooth (C°° or C w )-normal form is given by D = Ker a, 

y 2 . , „...,,. 

where a = dz dx. In this normal form we have the following identification : 

2 6 

• Martinet surface M : y = 0. 

d 

• The singular trajectories are the solution of Z = — restricted to y = 0. 
This case is called the Martinet case. 

Case 3. (Codimension 3). We take a point Xo € M PI S and we assume that 
the point qo is a regular point of M. The analysis of shows that in this 
case we have two different C°°- reductions to a C^-normal form depending both 
upon a modulus m. The two cases are : 

1. Hyperbolic case. D = Ker a, a — dy + (xy + x 2 z + mx 3 z 2 )dz. In this 
representation the Martinet surface is given by : 

y + 2xz + 3mx 2 z 2 = 

and the singular flow in M is represented in the (x, z) coordinates by : 

x = 2x + (6m — l)x 2 z — 2mx 3 z 2 
i = — (2z + 6mxz 2 ) 
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We observe that is a resonant saddle and the parameter m is an obstruc- 
tion to the C°° -linearization. 

2. Elliptic case. D = Ker a, a — dy + (xy + — + xz 2 + mx 3 z 2 )dz. The 
Martinet surface is here identified to : 

y + x 2 + z 2 + 2>mx 2 z 2 = 
in which the singular flow is given by : 

a; = 2z — h 6mx z — 2?ra z 

o 

z = — (2x + 6mxz 2 ) 

Hence is a center and still m is an obstruction to C°°-linearization. The 
analysis of 44 shows that the singularity is C°-equivalent to a focus. 

We call the case 3 the tangential situation because D is tangent to the 
Martinet surface at 0. We must stress that it is not a simple singularity and 
moreover there are numerous analytic moduli. 



2.5 Accessibility set near a singular trajectory 

The objective of this Section is to recall briefly the results of ^2] which de- 
scribe geometrically the accessibility set near a given singular trajectory satisfy- 
ing generic assumptions (see also |42l I41p . 



2.5.1 Basic assumptions and definitions 

We consider a smooth single input smooth affine control system : 

q(t) = Fo(q(t))+u(t)F 1 (q(t)), q e U. 

Let 7 be a reference singular trajectory corresponding to a control u S L°°[0, T] 
and starting at t = Q from 7(0) = qo and we denote by (7,p,u), where p is an 
adjoint vector for the associated solution of the equations from the maximum 
principle. We assume the following : 

(HO) (7,^) is contained in the set Q = {z = (q,p) ; {{H , Hi}, H{\(z) ^ 0}, 7 
is contained in the set ST = {q ; X(q) and Y(q) are linearly independant} 
and moreover 7 : [0, T] — > U is one-to-one. 

Then according to the results of Section 2.2, the curve z — (j,p) is a singular 
curve of order 2 solution of the Hamiltonian vector field H, with H = Hq + 

Yf — — — -j-^ — -j- Hi. Moreover the trajectory 7 : [0, T] — ► fi' is a smooth curve 
{{Hi, Ho}, Hi} 

and 7 is a one-to-one immersion. 

Using the feedback invariance of the singular trajectories we may assume 
the following normalizations : u(t) = for t E [0, T] and 7 can be taken as 
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the trajectory : t — > (i, 0, 0, ... ,0). Since w is normalized to 0, by successive 
derivations of the constraints Hi — 0, i.e < p(t), Fi(7(i)) >= for t <E [0, T], 
we get the relations : 

<p(i),V k (~f(t)) >=0 , fc = 0,...,+oo 

where T^ fc is the vector field ad k Fo(F\) and ad fe is defined recursively by : 

ad° F (Fi) = F u ad k F (Fi) = [F , ad^^o^i)] 

It is well known, see [2T], that for < > the space E{t) = Span jV fc (7(i)), fc G 
IN} is contained in the first order Pontryagin's cone K(t) evaluated along 7. We 
make the following assumptions : 

(HI) For t <E [0, t], the vector space i£(t) is of codimension one and generated 
b y {F°( 7 W),...,^"- 2 )(7(i))}. 

(H2) If n ^ 3, for each i e [0, T], X( 7 (t)) ^ Span {V°(7(i)), . . . , ^ ( "- 3) (7(t))} 

Definition 2.7. Let (^ (t) , p(t) , u(t)) be the reference trajectory defined on [0, T] 
and assume that the previous assumptions (HO), (HI), (H2) are satisfied. Ac- 
cording to (HI) the adjoint vector p is unique up to a scalar. The Hamiltonian is 
H = Hq + uHi along the reference trajectory and Hi = 0. If H = 0, we say that 

d dH 

7 is G-exceptional. Let D = ^^7^ =< P(t),[[Fi,F ],Fi](^(t)) > . The 
trajectory 7 is said G-hyperbolic if H.D > along 7 and G-elliptic if H.D < 
along 7. 

Remark 2.1. According to the higher- order maximum principle the condition 
H.D\ 7 called the Legendre-Clebsch condition is a time optimality necessary 
condition, see |23| . 

2.5.2 Semi-normal forms 

The main tool to evaluate the end-point mapping is to construct semi-normal 
forms along the reference trajectory 7 using the assumptions (HO, HI, H2) and 
the action of the feedback group localized near 7. They are given in |12j and we 
must distinguish two cases. 

Proposition 2.6. Assume that 7 is a G-hyperbolic or elliptic trajectory. Then 
the system is feedback equivalent in a C° -neighborhood 0/7 to a system (Nq,Ni) 
with : 



where a nn is strictly positive (resp. negative) on [0,T] if 7 is elliptic (resp. 

n—l q 

hyperbolic) and R = Ri-j—r is a vector field such that the weight of Ri has 



, <>''' 



order greater or equal to 2 (resp. 3) for i — 2, ... ,n — 1 (resp. i = 1), the 
weights of the variables q l being for i = 1, and 1 for i — 2, . . . , n. 
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Geometric interpretation : 

• The reference trajectory 7 is identified to 1 1 — ► (t, 0, . . . , 0) and the asso- 

d 



ciated control is u = 0. In particular N \~ = — — r 

1 ag! 



1 7 



• We have : 



d 



if k= l,...,n-2 



(ii) ad fc A^ .iVi| 7 = < 9g"- fe 

if fe > n - 2 

Ci") ad 2 7V 1 .AT = ^ 
and the first order Pongryagin's cone along 7 is : = 



2 ' ' 



i| 7 y| 7 

The linearized system is autonomous and in the Brunovsky canonical 
form : if 1 = ip 2 , . . . ,ip n = u. 

The adjoint vector associated to 7 is p — (e, 0, . . . , 0) where e = +1 in the 
elliptic case and e = — 1 in the hyperbolic case, the Hamiltonian being e. 

The intrinsic second-order derivative of the end-point mapping is identi- 
fied along 7 to : 



pi « 



aij(t) ip\t)ip>{t)dt 

with <^j 2 = </? 3 , . . . ,</? n_1 = <£™, y) n = u and the boundary conditions at 
s = and T : ip 2 (s) = ■ ■ ■ = (p n (s) = 0. 

Proposition 2.7. Let 7 be a G '-exceptional trajectory. Then n 3 and i/iere 
exists a C° -neighborhood of 7 m which the system is feedback equivalent to a 
system (N ,Ni) with : 



*=1 i,i=2 y 

« - 9 



Q q n-1 

" d 

where a n -i n -i is strictly positive on [0, Tl and R = > Ri^r^ , Rn-1 = is a 

t—f aq l 
i—i 

vector field such that the weight of Ri has order greater or equal to 2 (resp. 3) 
for i = 1, . . . , n — 2 (resp. i = n), the weights of the variables q l being zero for 
i = 1, one for i = 2, . . . , n — 1 and two for q n . 
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Geometric interpretation 



• The reference trajectory 7 is identified to t 1 — ► (t, 0, . . . , 0) and the asso- 
ciated control is u = 0. 



• We have the following normalizations : 

d 



(i) &d k N .N lh = \ dq?- 1 - 1 

for k > n - 2 



for k = 0, ... ,n — 3 



d 

(ii) A^o| 7 = ad n ~ 2 N .N lh = 



h 



a2 at 

(hi) ad 2 JV 1 .iVo= ^— ^ 



r ^ 

and the hrst order Pontryagin's cone along 7 is K\~ = Span < 
in the exceptional case j(t) is tangent to K\^uy 

The linearized system along 7 is the system : ip 1 — ip 2 , . . . , tp n ~ 2 = ip 
tp n ~ x = u. 



n-l 



„n-l 



• The adjoint vector p associated to 7 can be normalized to p = (0, . . . , 0, — 1). 

• The intrinsic second-order derivative of the end-point mapping along 7 is 
identified to : 

r T n—i 



/ a ^ (*) ^(*) ^'(*) * 

- 70 i,i=2 



with : ip 1 = (p 2 , . . . ,ip n 2 = ip n 1 , ip n = u and the boundary conditions 
at s = and T : <^(s) = • • • = p n-1 (» = 0. 

2.5.3 Evaluation of the accessibility set near 7 

We consider all trajectories q{t, u) of the system starting at time t = from 
7(0) = ; the accessibility set at time t is the set : A(0,t) = q{t,u). It is 

ueu 

the image of the end-point mapping. 

We use our semi-normal forms to evaluate the accessibility set for all trajec- 
tories of the system contained in a C°-neighborhood of 7. We have the following, 
see ^21 for the details. 
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Hyperbolic-elliptic situation 

By truncating the semi-normal form and replacing q 1 by t we get a linear- 
quadratic model : 

n 

q 1 = 1+ 5^ay(t)?V 
q 2 = q 3 , ... ,q n =u 

and it can be integrated in cascade. 

Let < t ^ T and fix the following boundary conditions : q(0) = and : 
<7 2 (t) = • • • = q n_1 (i) = = 0, we get a projection of the accessibility set 

A(0, t) in the line q 1 which describes the singularity of the end-point mapping 
evaluated on u(s) = for < s ^ t. Fig. represents this projection when t 
varies. 




elliptic case hyperbolic case 



Figure 1: 



Geometric interpretation The reference trajectory 7 is C°-time minimal 
(resp. time maximal) in the hyperbolic case (resp. elliptic case) up to a time t\ c 
which corresponds to a first conjugate time t\ c > along 7 for the time minimal 
(resp. time maximal) control problem. 

In particular we get the following Proposition : 

Proposition 2.8. Assume T < t\ c . Then the reference singular trajectory 7 
defined on [0,T] is in the hyperbolic (resp. elliptic) case the only trajectory 7 
contained in a C° -neighborhood 0/7 and satisfying the boundary conditions : 
7(0) = 7(0), j(T) — j(T) in a time T ^ T (resp. T ^ T). 
This property is called C° -one-side rigidity, compare with 
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Exceptional case 

We proceed as before. The model is 



q 1 = 1 + q 2 , q 2 = q 3 , . . . ,q n 1 = u 
n-1 

i,j=2 

Let < t,t' ^ T and consider the following boundary conditions : q(Q) = 
and = <, 9 2 (<') = ■ ■ ■ = q n_1 (t') = 0. We get a projection of the accessibil- 

ity set A(0,t') on the line q n . It is represented on Fig. [3] 




Figure 2: 



Geometric interpretation The reference trajectory 7 is C°-time optimal up 
to a time t\ cc which corresponds to a hrst conjugate time t\ cc > 0. In particular 
we have the following result. 

Proposition 2.9. Assume T < £ lcc . Then the reference singular exceptional 
trajectory 7 is C° -isolated (or C° -rigid). 

2.5.4 Conclusion : the importance of singular trajectories in optimal 
control 

The previous analysis shows that singular trajectories play generically an im- 

portant role in any optimal control problem : Min / f (x,u)dt when the 

Jo 

transfert time T is fixed. Indeed they are locally the only trajectories satisfying 
the boundary conditions and hence are optimal. If the transfert time T is not 
fixed only exceptional singular trajectories play a role. In fact as observed by [S] 
they correspond to the singularities of the time extended end-point mapping : 
E : (T,u) 1 — ► q(T,XQ,u). It is the situation encountered in sub-Riemannian 
geometry. 
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3 Generalities about sub-Riemannian geometry 



From now on, we work in the C w -category. 

Definition 3.1. A SR-manifold is defined as a n-dimensional manifold M to- 
gether with a distribution D of constant rank m n and a Ricmannian metric 
g on D. An admissible curve t i — ► q(t), ^ t ^ T is an absolutely continuous 
curve such that 4(0 £ ■^Gz(*))\{0} f° r almost every i. The length and the 
energy of q are respectively defined by : 

L{q)= f T (ffl,m) 1/2 dt, E{q)= f \q{t),q{t))dt 
Jo Jo 

where ( , ) is the scalar product defined by g on D. The SR-distance between 
qo,qi £ M denoted dsn(qo,qi) is the infimum of the lengths of the admissible 
curves joining q to q\. 



3.1 Optimal control formulation 

The problem can be locally restated as follows. Let qo £ M and choose a 
coordinate system (U, q) centered at qo such that there exist m (smooth) vector 
fields {Fi, . . . , F m } which form an orthonormal basis of D. Then each admissible 
curve t — > q(t) on U is solution of the control system : 

m 

q(t) = Y,<W(<i(t)) (ii) 

1=1 

The length of a curve does not depend on its parametrization, hence every 
admissible curve can be reparametrized into a lipschitzian curve s i — ► q(s) 
parametrized by arc-length : (q(s),q(s)) = 1, see |29| . 

If an admissible curve on U : t i — > q{t),0 ^ t ^ T is parametrized by 
arc-length we have almost everywhere : 

m tci 

q(t)=Y. U * F ^)> E U ?W = 1 

i=l i=l 

and L(q) = J Q T (J2 u i) 1 ^ 2 dt = T. Hence the length minimization problem is 
equivalent to a time-optimal problem for system This problem is not convex 

m 

because of the constraints : / JJ% = 1, but it is well-known that the problem 

i=i 

is equivalent to a time optimal control problem with the convex constraints : 

m 
i=l 

It is also well-known that if every curve is parametrized on a fixed interval 
[0, T], the length minimization problem is equivalent to the energy minimization 
problem. 
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Introducing the extended control system : 

m 

q(t) = ^u^FMt)) 

1=1 

m 

q°(t) = J2 u i®> ?°(0) = (12) 

and the end-point mapping E of the extended system : u ElA i — ► q(t, u, qo) , 
q = (q,q°),q(0) = (qo,0), from the maximum principle the minimizers can be 
selected among the solutions of the maximum principle : 

x, . dH x,, . dH dH 

*W = W p(t) = - W *r =0 

m m 

where H =< p, ^ \j,jFi(q) > +po^^uf is the pseudo-Hamiltonian and p — 

i=l i=l 

(p,Po) & ]R n+1 \{0} is the adjoint vector of the extended system. From the 
previous equation t i — ► po(t) is a constant which can be normalized to or 

m m 

— 1/2. Introduce =< p, Y^UjFj(g) > where ^ = or —1/2 ; then 

i=l i=l 

the previous equations are equivalent to : 

. _ dHu - = _dH^ dH v _ Q 
dp ' dq 1 du 

The solutions of these equations correspond to the singularities of the end- 
point mapping of the extended system and are called geodesies in the framework 
of SR-geometry. 

They split into two categories according to the following definition. 

Definition 3.2. A geodesic is said to be abnormal if v = 0, and normal if 
v = —1/2. Abnormal geodesies are precisely the singular trajectories of the 
original system (1 1 1 1) - 

A geodesic is called strict if the extended adjoint vector (p,po = v) is unique 
up to a scalar, that is corresponds to a singularity of codimension one of the 
extended end-point mapping. 

3.2 Computations of the geodesies 
3.2.1 Abnormal case 

They correspond to v = 0, and are the singular trajectories of system (|llfl . The 

m 

system is symmetric and hence H = V^UjP,, with Pj =< p, F. t (q) > . Therefore 

i=i 

they are exceptional. When m = 2, they are computed using the algorithm of 
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Section 2. The case m > 2 will be excluded in our forthcoming analysis because 
from [5] in order to be optimal a singular trajectory must satisfy the following 
conditions, known as Goh's conditions : 

<p(t),[Fv,Fw}(q(t))>=0 (14) 

Vu,u> S H m , Vf e [0,T] and Pm denotes ^UiFi. If m = 2, this reduces to the 
condition {Pi, P2] — deduced from the conditions P\ = P2 = but if m > 2 
it is a very restrictive condition which should not be generic (conjecture |21p. 

3.2.2 Normal case 

They correspond to v = —1/2. If the system of the Pi's is orthonormal then 







and hence Uj = P$ and P^ reduces to P„ = - ^~\P?, The trajectories 



au 2 

parametrized by arc length are on the level set H n = 1/2 and the normal 
geodesies are solutions of the following Hamiltonian differential equations : 

9= ^' P= -^a- (15) 

On the domain chart U, we can complete the m- vector fields {Pi, . . . , P m } to 
form a smooth frame {Pi, . . . ,F n } of TU. The SR-metric g can be extended 
into a Riemannian metric by taking the system of the Pi's as an orthonormal 
frame. We set Pj =< p, Fi(q) > for i = 1, . . . , n and let P = (Pi, • • • , P„). In 
the coordinate system (q, P) the normal geodesies are solutions of the following 
equations : 



m 

P = {P i ,Hn} = ^ l {P i ,P j }Pj (16) 

We observe that {Pi, Pj} =< p, [Pi, P/Kg) > and since the Pi's form a frame 
we can write : 



fi,fi](?)=E4(^(9) 



fc=i 



where the c,,'s are smooth functions. 



3.3 Exponential mapping - Conjugate and cut loci 

Assume that the curves are parametrized by arc-length. If t 1 — ► q(t) is any 
geodesic, the first point where q(.) ceases to be minimizing is called a cut-point 
and the set of such points when we consider all the geodesies with q(0) = qo 
will form the cut-locus L(qo). 
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The sub-Riemannian sphere with radius r > is the set S(qo,r) of points 
which are at SR-distance r from go- The wave front of length r is the set 
W(qo,r) of end-points of geodesies with length r starting from q . If Da.lXqo) 
is of rank n where Da.l is the Lie algebra generated by D, then according to 
Filippov's existence Theorem |27| if r is small enough each point of distance r 
from go is the end-point of a minimizing geodesic and S^go, r) is a subset of 
W(qo,r). We fix g £ U and let (q(t, qo,Po),p(t, qo,Po)) be the normal geodesic, 
solution of (|15|) and starting from (qo,po) at t = 0. The exponential mapping is 
the map : 

ex P 9o : (po,t)' — ►g(*,Po ) ?o) 

m 

Its domain is the set C x H where C is {po ; J^-P^Cpo, Qo) = 1}. If m < n it is 

i=i 

a (non compact) cylinder contrarily to the Riemannian case : m — n, where it 
is a sphere. 

A conjugate point along a normal geodesic is defined as follows. Let (po,ti) 
with t\ > be a point where exp go is not an immersion. Then t\ is called a 
conjugate time along the normal geodesic and the image is called a conjugate 
point. The conjugate locus C(go) is the set of first conjugate points. 

3.4 Gradated normal form 

3.4.1 Adapted and priviliged coordinate system 

Let (U, q) be a coordinate system centered at go, with D = Span {F±, . . . , F m }. 
Assume that D satisfies the rank condition on U. We define recursively : Do = 
{0},!) 1 = D and for p > 2 DP = Span {DP- 1 + [D 1 1 D p - 1 ]}. Hence DP is 
generated by Lie brackets of the Fi 's with length p. At p we have an increasing 
sequence of vector sub-spaces : {0} = D°(q) C D 1 (q) C • • • C D r ^ where r(q) 
is the smallest integer such that D r ^(q) = T q U. 

Definition 3.3. We say that go is a regular point if the integers n p (q) = 
dim D p (q) remain constant for q in some neighborhood of go- Otherwise we 
say that go is a singular point. Consider now a coordinate system (g 1 , . . . , g") 
such that dqi vanishes identically on D Wj ^ 1 (qo) and doesn't vanish identically 
on D Wj (go) for some integer Wj. Such a coordinate system is said to be adapted 
to the flag and the integer uij is the weight of q 3 . 

Definition 3.4. Consider now a SR-metric (D,g) defined on the chart (U, q) 
and represented locally by the orthonormal vector fields {F\, . . . , F m }. If / is a 
germ of smooth function at go, the order of / at go is : 

(i) if /(ffl,)^0, /*(/) = 0,ji(0) = +oo ; 

(ii) otherwise : fi(f) = inf {p / 3 Vx, ■ ■ ■ , V v G {Fi, . . . , F, n } with Ly t o ■■■ o 
Lv p (f)(qo) 7^ 0} where Ly denotes the Lie derivative. The germ / is 
called privileged if (J,(f) = min{p ; df(D p (qo)) ^ 0}. A coordinate system 
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{q l , . . . ,q n } is said to be privileged if all the coordinates qi are privileged 
at q . 

We have the following very important estimation, see [7] , [23] : 

Proposition 3.1. // (M,D,g) is a SR-manifold there exists a privileged coor- 
dinate system q at every point qo = of M. If Wi is the order (or weight) of the 
coordinate q d we have the following estimation for the SR- distance : 

d SR (0, (q\ . . . ,q n )) ~ l^l 17 ™ 1 + • • • + \q n \ 1 ^. 

Definition 3.5. Let (U, q) be a privileged coordinate system for the SR-structure 
given locally by the m-orhonormal vector fields : {F\, . . . ,F m }. If Wj is the 

, d 

weight of q J , the weight of — — is taken by convention as —Wj. Every vector 

oqt 

field Fi can be expanded into a Taylor series using the previous gradation and 
we denote by Fi the homogeneous term with lowest order —1. The polysystem 
. . . , F m } is called the principal part of the SR-structure. 

We have the following result, see 

Proposition 3.2. The vector fields F{,i = l,...,m generate a nilpotent Lie 
algebra which satisfies the rank condition. This Lie algebra is independant of 
the privileged coordinate system. 



3.4.2 Gauge classification 

Given a local SR-geometry (U, D, g) represented as the optimal control problem : 

m 

q = ^2 r u- l F l (q) 

i=l 
-T rn 

min / (j2uUt))dt, 
«(•) Jo K ~[ ' 

there exists a pseudo-group of transformations called the gauge group which is 
the subgroup of the feedback group defined by the following transformations : 

(i) germs of diffeomorphisms ip : q i — ► Q on U, preserving q ; 

(ii) feedback transformations u = (3(q)v preserving the metric g i.c, (3{q) € 
8(m,H) (orthogonal group). 

The invariants of the associated classification problem are the geodesies. 
They split into two categories : abnormal geodesies which are feedback invari- 
ants and normal geodesies. 

If q is an adapted coordinate system, a gradated normal form of order p ^ — 1 
is the polysystem {Ff , . . . , Ff n } obtained by truncating the vector fields F% at 
order p using the weight system defined by the adapted coordinates. 
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4 The role of abnormal minimizers in SR Mar- 
tinet geometry 



In this Section we analyze the role of abnormal minimizers in SR Martinet 
geometry which is the prototype of the generic rank 2 situation. Before to 
present this analysis it is important to make a short visit to the contact situation 
in H 3 . 



4.1 The contact situation in M 

The contact situation in H 3 has been analyzed in details in several articles 
|i3 0] . This analysis is based on computations about the exponential mapping 
using a gradated normal form. To understand the remaining of this article it is 
important to make the contact situation fit into the following framework. 

First, without losing any generality we can use to understand a generic 
contact SR-problem a gradated form of order 1 computed in 0] where the SR- 
metric is defined by the two orthonormal vector fields : F±, F 2 where : 

d y . d d x , n . d 

where Q is a quadratic form : ax 2 + 2bxy + cy 2 depending on 3 parameters. 
The weight of x, y is one and the weight of z is two. When a = b = c = 0, 
it corresponds to the contact situation of order -1 which is the well-known 
Heisenberg case but also a gradated normal form of order 0. 

d 

To get an adapted frame we complete F\ , F2 by F3 = — . Computing we 

dz 

d 

get : [Fi,F 2 ] = (1 + 2Q)— . Using Pj =< p,Fi(q) >, the geodesies equations 



are 



x = Pi 

y = P2 

Piy(l + Q) P 2 x(l + Q) 

Z = 2 2 

A = {Pi,P 2 }P2 = (1 + 2Q) P 2 P 3 

Pi = {P2,Pi}Pi - -(1 + 2Q)P 1 P 3 

P 3 = 



In the Heisenberg case we have Q = 0, and if we set P3 = A we get : 
Pi + A 2 Pi =0, which is a linear pendulum. 

Using the cylindric coordinates : P± = sin#, P 2 = sin^jPj = A, where 
9 ^ kir, the geodesies parametrized by arc-length are solutions of the following 



22 



equations : 

. . P lV {\ + Q) P 2 x(l + Q) 
x = Pi y = P 2 z= 

(17) 

9 = (1 + 2Q)A 

where A is a constant. The important behavior is when A — > oo. We may assume 
A > 0. By making the following reparametrization : 

ds = A(l + 2Q)dt (18) 
d9 

the angle equation takes the trivial form : — = 1. Hence it is integrable and 

we obtain 9(s) = s + 8 - 

The remaining equations take the form : 

dx sin 9{s) 



(19) 



ds (1 + 2Q)A 

dy cos 9(s) 

ds ~ (1 + 2Q)A 
dz sm9{s){y{\ + Q)) - cos 9(s)(x(l + Q)) 

~d~s ~ (1 + 2Q)A 

For large A, they can be integrated as follows. We set e = 1/A : small parameter, 

x = eX, y = eY, z = e 2 Z, - — \— = 1 + Q = 1 + Ax 2 + 2Bxy + Cy 2 H and 

1 

we get : 

X = sin (s + 9 Q )[l + e 2 Q(X,Y)+o(e 2 )} 
Y = cos (s + ) [1 + e 2 Q{X, Y) + o(e 2 )] 

• sin (s + 9 )Y - cos(s + 9 )X 
Z = +o(e) 

The previous equations can be integrated by quadratures by setting : 

X = X + e 2 X 1 +o(e 2 ) 

Y = Y a + e 2 Y x + o(e 2 ) 

Z = ^ + o(e) 

and we get in particular 

X a = sin(s + ) 
Y = cos(s + 9 ) 

Zo 



sin(s + e )Y (s) - cos (s + 9 ) X (s) 



2 

Xx = sin(s + o ) Q(X ,Y a ) 
Y = cos{s + 9 ) Q(X Q ,Y ). 
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The solutions are computed in the s-parametrization and the arc-length t can 
be computed by integrating (|19fl by quadratures. 

If we want to mimic this procedure in the Martinet situation, we shall en- 
counter integrability obstructions due to the existence of abnormal geodesies. 

The sphere in the flat contact situation is represented on Fig. [3J 




Figure 3: SR sphere in the flat contact case 



4.2 The Martinet situation 
4.2.1 Normal forms and invariants 

The Martinet SR-geometry is rather intricate and it is difficult to make a priori 
normalizations. It will appear later that a good starting point to make the 
computations is to use the following normal form computed in |2] : 

• The distribution D is taken in the Martinet-Zhitomirski normal form : 

, y 2 

D = Ker lu, lo = dz dx. 

2 

• The metric on D is taken as a sum of squares : a{q)dx 2 + c(q)dy 2 . 

In this representation the Martinet surface containing the abnormal geodesies 
is the plane : y — and the abnormal geodesies are the straight-lines : z = z Q . 
The abnormal line passing through is given by 7 : 1 1 — > (±£, 0, 0). 

The computations in [2] show that we can make an additional normalization 
on the metric by taking either the restriction of a or c to the Martinet plane 
y = equal to 0. 

The variables are gradated according to the following weights : the weight 
of x, y is one and the weight of z is three. By identifying by convention at order 
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p two normal forms where the Taylor series of a and c coincide at order p we 

end up with the following representatives of order : 

either 

g=(l + ayfdx 2 + (1 + /3x + -yy) 2 dy 2 

or 

g = (1 + ax + j3y) 2 dx 2 + (1 + jy) 2 dy 2 . 

In each of those representations the three parameters are, up to sign, invariants. 
They can be used to compute the exponential mapping in the generic situation. 
If we truncate g to dx 2 + dy 2 it corresponds to the principal part of order — 1 of 
the SR-structure defined previously In the sequel it will be called the flat case. 

4.2.2 Geodesies equations 

The distribution D is generated by : 

r< d y 2 d d 
Cxi = h — -7T and G 2 = — 

ox 2 oz oy 

and the metric is given by g = adx 2 + cdy 2 . We introduce the frame : 

*i = -U, F 2 ='G 2 , F 3 = A 
and Pi —< p, Fi(q) > for i = 1, 2, 3, i.e 

p _ Px+PzV 2 /2 Py 

First, we assume that g is not depending on z ; this is the case for the 
gradated normal form of order 0. It corresponds to an isoperimetric situation, 

d 

that is the existence of a vector field Z identified here to — transverse at to 

oz 

D(Q) and the metric g does not depend on z. 
The system is written : 

. y 2 

x = ui , y = u 2 , z = —ui 

and the Hamiltonian associated to normal geodesies is : 

H n {q,p) = ^(uja + ulc) 



and the geodesies controls are : 



til = - (Px+P z y 2 /2) , u 2 = — • 
a c 
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Normal geodesies are solutions of the following equations : 

x = -(px +Pz y 2 /2), y = — , z = ^-{px+Pz y 2 m 

a c la 

. _ (Px+Pz y 2 /2) 2 

Px 2c 2 + 2a 2 Gx (20) 

P 2 v c v . {Px+Pz y 2 /?) 2 {Px+Pz y 2 /2) 

Py = + 2^ ay a PzV 

Pz=0 

In the (q, P) representation the previous equations take the form : 

x = ^- i = — z= V —^ 
\fa \fc 2 yfa 

Pi ( Oril C X 



v^v^ V 2ja~ 2^fc V (21) 
■ P\ ( a„ c x 

P2 = —r±7= (yPz - ttfPi + irk ft 

ay/c V 2y/a zyc 



P 3 =0 

If we parametrize by arc-length and if we introduce the cylindric coordinates : 
Pi = cos 9, P2 — sin 9, P3 = A, we end up with the following equations : 

P_ 

x — , 

\ a 



1 



v a v c 

Pi = X 



P2 


■ - y2 Pl 




2 y/E 


yPs- 


dy 



c x 

P2 



2^fc 



(22) 



It is proved in that for a generic SR-problem, each geodesic is strict. In 
our representation we have the following result. 

Lemma 4.1. The abnormal geodesic 7 : t 1 — ► (±t, 0, 0) is strict if and only if 
the restriction of a y to the Martinet plane y = is 0. 

Using the gradated normal form of order with the normalizations : 
a=(l + ay) 2 , c = (1 + (3x + jy) 2 
the equations 11' 2 1) reduce to : 

cos 9 sin 9 y 2 cos 9 

x = —r^ y=~r z = ^ — t- 

V a V c 2 V a (23) 



[yX — acos 9 + /3sin 
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The previous equation defines a foliation (!F) of codimension one in the plane 

(y,ff). Indeed using the parametrization : yfa \fc — = — — and denoting ' the 

ere dr 

derivative with respect to r, the equations can be written : 
x = \fc cos 9 z = \fc — cos 9 

(24) 

y' = \J~a sin 9 9' = —[y\ — acos 9 + /3sin 9] 

and they can be projected onto the plane (y,9). The last two equations are 
equivalent to : 

9" + Asin 9 + a 2 sin 6>cos 6 - a(3sin 2 9 + /30'cos 9 = 0. (25) 

This equation will be used in the sequel to study the SR-Martinet geome- 
try in the generic case of order 0. Unfortunately it depends on the choice of 
coordinates. Note that in the flat case where a = c = 1 the equation reduces to 
9" + A sin 9 = which is a nonlinear pendulum. 

4.2.3 Conservative case 

The analysis of Subsection 4.1 shows that in the contact case the equation l|17fl 
associated to the evolution of 9 defines an intcgrablc foliation. In the Martinet 
case the foliation defined by equation (|25|l is not in general integrable. This leads 
to the following definition which is independant of the choice of coordinates. 

Definition 4.1. Let e(t, 9, A) be a normal geodesic parametrized by arc- length 
starting from q(0) = and associated to 9(0) = 9 , f*3(0) = A. The problem is 
said conservative if there exists a coordinate y transverse to the Martinet surface 
such that for a dense set of initial conditions (9 , A) the trajectory t — > y(t) is 
periodic up to reparametrization. The equation describing the evolution of y is 
called the characteristic equation. 

4.2.4 Analysis of the foliation T 

The foliation (F) is described by equation (|25|l : 

9" + Asin 9 + a 2 sin 9 cos 9 - af3 sin 2 6> + f3 9' cos 9 = 0. 

Moreover recall the relation : y' = (1 + ay)cos 6, 9' = — (y\ — acos 9 + (3 sin 9) 
The singular line project onto 9 = kit which correspond to the singularities 

of : 9 = kn. & = 0. ' 

Among the solutions of (|25|l . only those satisfying the relation : 

6' = a cos 9 + /3 sin 9 (26) 

at t = correspond to projections of geodesies starting at t = from q(0) = 0. 
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Using an energy-balance relation we can represent the solutions of (J 7 ) for 
| A \^>\ a |,| |,| 7 |, see |TJ]]- We may suppose A > 0. Introducing the small 
parameter : e = 1/ \/A and the parametrization s = tVX we get the equation : 



£0 

ds 2 



sin 



e(3cos 



d6 
ds 



e 2 a sin 9(a cos 9 - /3 sin 6») = (27) 



and equation l(2*B|) takes the form : 

e(a cos + (3 sin 



The flat case corresponds to a = /3 = 0, i.e 

s = and is aZso the limit case e — ► 0. 
The following result is straightforward. 



ds 2 



Sill 



(28) 



de 

= 0, — = at 
ds 



Lemma 4.2. The problem is conservative if and only if [3 = 0. 

We represent below the trajectories of (J 7 ) for A a |, | j3 |, | 7 |, on 
the phase space : (9, 9) but geometrically it corresponds to a foliation on the 
cylinder : (e %e , 9). 

• Flat case (a = j3 = 0). It corresponds to a pendulum, see Fig. 0] 




Figure 4: 

The main properties are the following. We have two singularities : 

— is a center. 

— (7r, 0) is a saddle and the separatrix S is a saddle connection. 
Only the oscillating trajectories correspond to geodesies starting from 0. 
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Conservative case ((3 — 0) The equation reduces to 
ds 2 



— — + sin 9 + e 2 a 2 sin 9 cos 9 = 



d9 

Multiplying both sides by — and integrating on [0, s] we get 

ds 



1 /ff(9\2 

2 \ds 



cos 9 



2 2 2 

e a cos 



and the system has a global C u first integral 

1*. 



V(9,9) = -9 2 - (cos 



■cos 2 (9 



(29) 



The phase portrait is similar to the one in the flat case but the section 

d9 

defined by (|2*5|l and corresponding to y = is here : — = eacos 9. In 

particular if a ^ (strict case) there exist both oscillating and rotating 
trajectories corresponding to projections of geodesies starting from 0, see 
Fig.0 



ds 




Figure 5: a > 

• General case (/3 ^ 0) The two main differences are the following : 

— the center becomes a focus ; 

— the saddle connection is broken. 

The trajectories are represented on Fig. [S] 

The respective generic behaviors of 1 1 — > y(t) are represented on Fig. [7| 

It is important to observe that our description of the behavior oft i — ► y(t) 
is true for the gradated form of order 0, but also of any order when X — > 
oo. 
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Figure 6: > 




Figure 7: ^ 
4.2.5 Characteristic equation 

If = 0, the Hamiltonian iJ„ = ^(-P 2 + ^2) nas two cyclic coordinates : x 
and z and therefore = cos 9(0) and = A are first integrals. The equation 

H n = 1/2, with Pi = Px + P f. y ^ and P 2 = % takes the form : 

V a v c 

2 /9x 2 



dt . , 

Introducing : dr = it becomes 



(f ) 2 + (^+P,2/ 2 /2) 2 = a 
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where a = (1 + ay) 2 . Hence we get 



iL) 2 = F{y) 00) 



At, 

where F(y) = (1 + ay) 2 — (p x + p z y 2 /2) 2 . The analysis is based on the roots of 
the quartic F{y). We assume A > 0. 

We observe that F can be factorized as F1F2 with : 

F x = (1 + ay) - (p x + p,y 2 /2) , F 2 = (1 + ay) - (p z + p 2 y 2 /2) 

and we can write : 

F{y) = (2m 2 -±(y-^) 2 ){2 m " + \{y + ^) 

where : 2m 2 = 1 — p x + ^ , 2to" = 1 + p x — |^ 

and : m 2 + m" = 1 , m 0. 

If we set : 77 = « ^ = + _« 

' 2m 2mVA ' ' 2m 1 2mVX 

we can write : 

F(y) = 4m 2 (l-r) 2 ){m" + m 2 f) 2 ) (31) 

\l to!* 

F is a quartic whose roots on C are r\ = ±1 , fj = ± . 

m 

The case m" = is called critical and it corresponds to a double root for F. 
We have : 

Lemma 4.3. In the strict case a ^= 0, there exist geodesies starting from 
which are critical. 



Geometric interpretation 

The critical geodesies project in the (6>, 9) phase space onto a separatrix, see 
Fig. 

The characteristic equation can be put into a normal form using an homo- 
graphic transformation to normalize the roots of F. The procedure is standard, 
see [2S]- We proceed as follows ; F is factorized into F1F2 and we consider the 
pencil F\ + vF 2 of two quadratic forms. If a. 0, there exist two distinct real 
numbers v\, v 2 such that F± + vF 2 is a perfect square : K\{y — p) 2 , K%(y — q) 2 . 
Using the homographic transformation : 

V-P /o \ 
u = , 32) 

y-q 

the characteristic equation can be written in the normal form : 

dy (p — q) _1 du 



/F(y) V(^ 2 + Bi)(A 2 u 2 + B 2 )' 



(33) 



The right hand side corresponds to an integrand of an elliptic integral of the 
first kind. More precisely, excepted the critical case m" — 0, the solution y in 
the w-coordinate can be computed as follows : 
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• if the quartic F admits two real roots, u can be parametrized using the 
cn Jacobi function ; 

• if the quartic F admits four real roots, u can be parametrized using the 
dn Jacobi function. 

If a = 0, the analysis is simpler, indeed F(y) can be written : 
F(y) = 4fc 2 (l - 77 2 ) {k' 2 + k 2 n 2 ) 

where 77 = ^/—JL and 77 can be computed using only the cn function. 
2k 

Proposition 4.4. We have two cases : 
2k 

(i) If a — 0, y — —/=n where 77 is the cn Jacobi function. 

V A 

(ii) If a =/= 0, y is generically the image by an homography of the cn or dn 
Jacobi function. 

Geometric interpretation If a = 0, the motion of y is a cn whose amplitude 
2k 

is —= . The motion is symmetric with respect to y = and the amplitude tends 
vA 

to when A tends to the infinity, see Fig. |8] 

ua — p 

If a =/= 0, we can expand : y — near u — 0. The motion of y is no 

u — 1 

more symmetric with respect to y — and there is a shift. Hence y can be 
approximated by a constant plus a cn or dn motion. 




Figure 8: a > 
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4.2.6 Integral formulas in the general conservative case 

If the metric g does not depend on a;, it is convenient to use the following integral 
formulas from |24j to compute x and z in terms of y. 

We denote by e(t),t £ [0, T] a normal geodesic starting from and we 
assume that the component : t i — ► y(t) oscillates periodically with period V. 
We denote by < £i < ■ ■ ■ < ^ T the successive times such that y{U) = 0. 
We introduce : 

sign y(0) if 2/(0) + 
sign 2/(0) if 2/(0) = 



and we set : 



V+= max 2/(i) , y- = min y(i) 
te[o,P] te[o,F] 



Parametrizing the geodesies by y we must integrate the equations : 



/c 



dx dz_ yl^/l^L dt- - // 

dy ~ VoP a ' dy ~ 2 v^P 2 ' ~ P 2 y 



where P 2 (y) = ay/l-P?(y) for i € [0, t{\. 

This allows to get explicit integral formulas. In particular if y(T) = for 
T = ijy we get : 

• N odd 

x(t) = 2 r ^ - dy +{ N-D r 4 , PM * 

y o v/i-^ 2 ^) y»- ^ 

v ° - r c y 2 Pi(y) J , „ T + y 2 A(y) 



r(T) = / a 
Jo 



_ dy + (N - 1) / a ya ' dy 

(34) 



• TV even 



*(T) =N T ^ £M dy 
Jy- 2v^ ^-^(y) 



(35) 



and the period is given by : 

V = 2 f J+ - ^ dy. (36) 

The integrands have simple poles when P\{y) = ±1 so the integrals are 
well-defined. 



33 



4.2.7 The return mapping 

The main geometric object to understand the role of abnormal trajectories in 
the problem is the return mapping. Indeed if we consider the trace of the sphere 
and the wave front in the plane y = : 

S(0, r) = 5(0, r) n (y = 0) , W(0, r) = W(0, r) n (y = 0) , 

they are in the image of the following mappings. 

Definition 4.2. Let e : (i G [0,T], 0(0), A) i — ► (x(t), y(t), z(t)) be a normal 
geodesic, parametrized by arc-length. If y{t) ^ 0, we can define < t\ < ■ ■ ■ < 
tN $J T as the times corresponding to y(U) = 0. The first return mapping is : 

Ri ■■ (A, 0(0)) G£>r^ (xfa^zfa)) 
and more generally the n-th return mapping is the map : 

R n : (A, 0(0)) G D n i — > (a;(i n ), z(i n )) 
where Z)j are the domains. 

If the length is fixed to r, we observe that W(0, r) is the union of the image 
of the return mappings and (±r, 0) which are the end-points of the abnormal 
geodesies. 

The remaining of this Section is devoted to the analysis of the return map- 
ping. We proceed by perturbations of the flat case. We shall estimate the 
asymptotic expansions of S and W in the abnormal direction. They are an 
union of curves in the plane. Such a curve is subanalytic if and only if it admits 
a Puiseux expansion. It is a practical criterion to measure the transcendence of 
the sphere and wave front in the abnormal direction. 

4.2.8 The pendulum and the elastica in the flat case 
In the flat case the equation 1)27(1 is a simple pendulum : 

^+sin0 = O 
ds z 

d9 

where s = t^/X, t is the arc-length parameter and y = -j= — . In particular if 

V Xds 

d9 

y(0) = 0, we have — = 0. We get : 
ds 

(_J = 2(cos - cos 0{0)) 

The integration is standard using elliptic integrals The characteristic 
equation takes the form : 

is 1 = (l-Px- y Pz ) i 1 + Px + ^2 Pz ) 
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and we introduce k, k' £ [0, 1] by setting : 

2fc 2 = l-p x , 2fc' 2 = l+p x 

where p x — cos 9(0). Wc set : r} = and we get the equation : 

^ = (1-ry 2 ) (fc' 2 + fcV) 

We integrate with rj(0) = y(0) = and we choose the branch ?)(0) > corres- 
ponding to y(0) = sin 9(0) > 0. We get using the cn Jacobi function : 

r/(t) = -cn (K(k)+tV\,k) 

where AK(k) is the period, K being the complete elliptic integral of the first 
kind : 

K(k) = f dv = (1 - k 2 sm 2 9)-U9 

Jo yJii-rfW 2 + k 2 r] 2 ) Ja 

Hence 

2k 

y(t) = -=ua.(u,k), u = K + tV\ (37) 

V A 

which coincides with the formula obtained by integrating the pendulum. 
The components y and z can be computed by quadratures and we get : 

x(t) = -t+^= (E(u)-E) 

V (38) 
z(t) = ^jj^ [(2k 2 - 1) (E(u) - E(K)) + k' 2 tV\ + 2fc 2 sn u cn u dn u] 

where E is the complete elliptic integral of the second kind : 

f K(k) ,77/2 

E(k) = dn 2 u du = (l-k 2 sin 2 9) 1 ' 2 d9 

Ja Jo 

The previous parametrization corresponds to geodesies with A > 0, 9(0) € 
]0, 7r[. The solutions corresponding to A > 0, 9(0) e] — it, 0[ are deduced using the 
symmetry : Si : (x,y,z) i — > (x,—y,z). The solutions corresponding to A < 
are deduced using the symmetry : 52 : (x, y, z) i — > (—x, y, —z). The solutions 
with A = play no role in our analysis. 
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•V 



•y 




(i) (ii) 
Figure 9: elastica 




Figure 10: behaviour on the separatrix 

Elastica 

The projections of the geodesies on the plane (x, y) are parametrized by : 

2k 2 
y(t) = -^cn (u, k) , x(i) = -t+ ^( E ( u ) - E ) 

They are precisely the inflexional elastica described in |31| . 

They take various shapes whose typical ones are represented on Fig. |5J 
When k! — > the limit behavior is represented on Fig. (ii) , see also Fig. ^| 
(behaviour on the separatrix). 

In this representation 9 is up to a constant the angle of the normal with 
respect to a fixed direction. The rotating trajectories of the pendulum corres- 
pond to geodesies not starting from 0. They project on the space (x, y) onto 
non inflexional elastica, see Fig. 1111 (ii). 
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(i) 



(ii) 



Figure 11: 



4.2.9 Trace of 5(0, r) and W(0,r) in y = in the flat case 

2K 

The successive intersection times with y = are given by : ti — —= , i = 1, 

V A 

If we fix the length to t, — r, we get the following curves : 



, AT. 



x = -r + —= (E(K + i2K) - E) 
V A 



3A 3 /2 



(2fc 2 - 1) (J5(i? + i2i^) - E) + 



It represents a parametric curve, where the parameter is k g]0, 1[. Using the 
relation : + i2K) = (2i + 1)E we obtain for each i the following curves : 

k i — > Cj(fc) = (xi(k),Zi(k)), 



Xi(k) 
Zi{k) 



E 

-r + 2r- 



6i 2 K 3 



(2k 2 - l)E + k' K 



where k €]0, 1[. We can easily draw those curves using the standard package 
about elliptic functions in Mathematica, see Fig. ^| 




(-r,0) 



(r,0) 



Figure 12: 
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The exterior curve obtained for i = 1 represents the intersection of the sphere 
5(0, r) with the Martinet plane in the domain z > 0. Each point of this curve 
is the end-point of two distinct minimizers and by obvious geometric reasoning 
we have : 

Proposition 4.5. The cut locus L(0,r) is C\ U —C\. 

Moreover by inspecting Fig. ^| we deduce the following : 
Proposition 4.6. The abnormal geodesies are minimizers. 

This result is not new but here the proof is based on the analysis of the 
geodesic flow. The main property is that at each intersection with y = 0, the 
variable z has non zero drift which can be easily evaluated using This will 

lead to optimality results for the general metric, by stability. 

This is an alternative proof to the optimality results presented in Section 3 
or in U2|, [221 , where we consider all the trajectories of the system. 

Remark 4.1. We observe that (— r, 0) is a ramified point of the trace of the wave 
front on the Martinet plane with an infinite number of branches. This gives us a 
precise geometric interpretation on the structure of the geodesies of fixed length 
with respect to the abnormal line. Indeed for every neighborhood U of (— r, 0, 0) 
and every n G M, there exists a geodesic of length r with end-point in U , with 
n oscillations. 

We represent on Fig. 1131 the first and second return mapping, the length 
being fixed to r, and by restricting the domain to A > 0, 0(0) G [0, n]. 




Figure 13: 

In the phase space (0,0), Ri corresponds to the symmetry : (0,0) i — ► (~ #,0) 
and i?2 corresponds to the identity : (0, 0) i — ► (0, 0). 

We represent on Fig. ^] the two branches C\ and G\ in 5(0, r) ending 
at (— r, 0) and (r, 0) and corresponding respectively to the behaviors of the 
geodesies near the center and the separatrix E. 

Inspection of Fig. ^] leads to the following. 
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Figure 14: 



Proposition 4.7. For each n ^ 1, the return mapping R n is not proper. 

Proof. The inverse image of a compact ball centered at (— r, 0) corresponds to 
an asymptotic branch in the parameter space (0(0), A). The transcendence of 
this branch can be easily computed. Indeed when k' — ► 0, K(k') ~ In 1/k' and 
the branch is logarithmic. □ 

4.2.10 Asymptotics of the sphere and wave front near (r,0) and (-r,0) 

We can estimate the branches C\ and C\. The computations are geometrically 
different. Indeed the computation of C\ requires the estimation of the leaves 
of the foliation J 7 , localized near the center but the computation of C\ requires 
the estimation of the leaves near the separatrix E connecting the saddle points 
(—it, 0) and (tt, 0). To make the estimation we use the parametric representation : 

x(k) = — r + 2r — 
K 

z(k)^^[(2k 2 -l)E + k' 2 K] 

where k g]0, 1[ and C% (resp. C\) is obtained by making k — > (resp. k — > 1). 

The transcendence of the branches is related to the properties of the complete 
integrals : 

K= f dV = C'\l-k 2 sin^)- 1 / 2 old 

Jo y / (i- r? 2)(fc'2 + fc ,2 + fc2?72 ^ Jo 



and 



I dn 2 udu= I (1 - k 2 sin 2 9) 1/2 d9 
Jo Jo 



Both E and K are solutions of hypergeometric equations whose singular 
points are located at k — and 1. Using this properties we deduce the following 
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Lemma 4.8. When k — ■+ 0, E and K are given by the following converging 
asymptotic expansions : 

Lemma 4.9. When k! = yl — k 2 — > we have : 

E(k) = Ul (k' 2 )ln 1 + u 2 (k' 2 ) 

K{k) =u 3 (k ,2 )ln - + u A {k' 2 ) 
k 

where the Ui 's are analytic near and can be written as : 

Ul (k' 2 ) = %- + 0(fc' 4 ) , u 2 {k' 2 ) = \- k — + 0(fc' 4 ) 

2 2 

u 3 (k' 2 ) = 1 + ^ + 0(fc' 4 ) , u 4 (fc' 2 ) = ~+ 0(k /4 ). 

Remark 4.2. The complete expansions are given in |18| . The general theory 
about Fuchsian differential equations guarantees the convergence of the previous 
expansions and the coefficients can be recursively computed using the ODE. 
Another method which can be applied in the general conservative case is to use 
the integral formulas. 



Estimation of C\ 

When k — > 0, E and l/K are analytic and we have the following estimations 
using Lemma f4. 81 : 

E k 2 

— = 1 +o(k 2 ) 

K 2 y ' 

l2 , „//„2\ 



x(k) - r = —rk + o(k ) 



z{k) = —k 2 +o{k 2 ) 



In particular we deduce the following : 

Proposition 4.10. When k — > 0, the branch C\ is semi-analytic and is given 
by a graph of the form : 

2r 2 

z = —-—^ (x — r) + o(x — r) , x ^ r 
6ir z 



Estimation of C\ 

When k! — > 0, we cannot work in the analytic category but in the log-exp 
category introduced in |19j . Using |28j . the elimination of the parameter k' is 
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allowed in this category and will lead to a log-exp graph. The precise algorithm 
to evaluate C\ has been established in [5] and we proceed as follows. 

X ~\~ V Z 

We set X = — — , Z = — , and we get : 
2r r A 

X=§, Z=^[(2P-l)E + k' 2 K] 

If we introduce : X\ = k! , X2 = ln \j k , , we have Xi , X2 — > when fc' — > + 
and both X and Z are analytic functions of X\ and X 2 . 
An easy computation shows that : 

Xi~4e - ^, X 2 ~X when 1^0+ 

and we can write : 

Xi = 4e~y (l + Yi(X)) , X 2 = X(1 + Y 2 (X)) 

where F x , r 2 -> when X -> 0+. 

Both Yi and Y2 can be compared and a computation gives us : 

Y 2 = XA ± (X, Y0 , Yi ~ ^ when X -> 0+ 

where is a germ of an analytic function at 0. 

Now the equation X = E/K can be solved in the variables Yi, Xi, X 2 using 
the Implicit Function Theorem in the analytic category and the computations 
show that : 

1 

Yl =A 2 {X, e -^-) 

where A 2 is a germ of an analytic function at 0. Using this relation we end 
with : 

1 

e x 

Z = F(X,—) 

where F is a germ of an analytic function at 0. 

This is the constructive algorithm to compute the branch C\ as a graph in 
the log-exp category. Hence Z can be expanded as : 



+00 



u p (X) 



^ 1 x 



To ensure that G\ is not semi-analytic we must check that there exists a non 
zero term of the form ak, p X k (e~^) p , p > in the expansion. For this we 
compute the first non zero coefficient according to the lexicographic order on 
the pair (p, k). The simplest computation made in |2] is to observe that : 



E .„ l-2fc' 2 / E\ 3 fc' 2 



X = — , 6Z 



K ' E 2 \K K 
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but the algorithm which can be generalized is the following. We use the approxi- 
mations : 

fc' 2 fc' 2 
E = 1 + — In 4/fc' - — + o(fc' 2 ) 
2 4 

fc' 2 

K = In 4/fc' +— In 4/fc' + o(fc' 2 ln 4/fc'). 
Easy computations lead to the formula : 

3 5 fc' 2 / fc' 2 \ 

~4 (In 4/fc') 4 + °Un4/fc'J 

Using fc' ~ 4e - ^, — \- rn ~ X we obtain : 

& 1 In 4//r 

Z = -X 3 - 4c~y x 3 + o(X 3 e"*) 
6 

Remark 4.3. We observe the following : 

• Uq(X) — X 3 /6 is algebraic. 

• There is a phenomenon of compensation and the first non zero flat term 
is of the form X 3 e~^ and not X 2 e~^ ; that's why we need three terms 
in E and two terms in K. 

• In general the computation of the first non zero a p ^ can be done in a 
finite number of steps, for instance using a finite number of coefficients of 
u (X). 

4.2.11 Numerical aspects 

Fig. E| represents the numerical simulation of the flat Martinet sphere. We 
observe a numerical problem when computing near the abnormal direction. 

4.2.12 Asymptotics of the sphere and wave front in the abnormal 
direction in the conservative case 

Geometric preliminaries 

We can estimate the sphere and the wave front in the abnormal direction 
when g = (1 + ay) 2 dx 2 + (1 4- 'yy) 2 dy 2 (or in the general case) using the integral 
formulas l|34l) . We observe that the geometry remains invariant for the following 
symmetry : S\ : (x, y, z) i — > (— x, y, —z) and in our study we can assume A ^ 0. 
Another symmetry is the following. Adding to the geodesies the equations : 
a = 0, 7 = we can observe that the geodesies equations are left invariant by 
the transformation : (x, y, z,p x ,p y ,p z , a, 7) 1 — > (x, -y, z,p x , -p y ,p z , -a, -7). 
Hence we can fix the sign of one of the parameters and we shall make the 
following choice : a ^ 0. 
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Figure 15: Flat Martinet sphere 



Let e(t) = (x(t), y(t), z(t)) be a normal geodesic starting from and asso- 
ciated to p y (0) = sin 0(0), p x = cos 0(0) and p z = A. We observe the following. 
If A is non zero the y component of a geodesic oscillates periodically unless 
it corresponds to a separatrix E between two values y- and y+ and we have 
y- < < y+ if y(0) 7^ 0. If y(0) = 0, then sign y(0) = sign a > when a > 0. 

Moreover using Fig. 0or the integral formulas 135(1 . we deduce the following 
Proposition. 

Proposition 4.11. Let e(t) = (x(t), y(t), z(t)) be a geodesic starting from 
such that y oscillates periodically, y(0) ^ and corresponding to the initial 
conditions y(0),p x and p z . Lete(t) = (x(t),y(t),z(t)) be the geodesic associated 
to — 2/(0), p x andp z . Then e ande are distinct but their even intersections with 
the plane y = are identical and have the same length. In particular e(.) is not 
a minimizer beyond its second intersection with the plane y = 0. 

This is illustrated on Fig. 1 161 where we project a geodesic in the plane (a;, y). 
Conclusion 

The previous Proposition tells us that except when p y (0) = sin 9(0) = 0, the 
sphere is contained in the image of i?i and R2 ■ The others cases can by studied 
by continuity or using a numerical algorithm developped in |17j to compute the 
conjugate points. 

We shall now estimate the image of R\ and R2 near the two singularities of 
the foliation T . 
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Figure 16: 

Estimation of R\ 

dO 

The constraint y = takes the form S : — = eacos where cos can 

be approximated by ±1 near 9 = 0,it. Contrarily to the flat case we must 
distinguish the case 9(0) e] — ir, 0[ where a = sign y(0) = +1 from the case 
0(0) e]0,7r[ where a = —1. We use following notations : 

• C(D) branches corresponding to an oscillating (resp. rotating) pendulum 
or CD : mixed behaviors. 

• Symbols without bars : behavior near the separatrix, symbols with bar : 
behaviours near the focus. 

• When a = +1, we use the symbol '. 

They are images by R\ of curves in the parameters A, 6(0) denoted by the same 
but minuscule symbol. We obtain the Fig. IT7I 

Estimation of R2 

The analysis is simpler because the branches corresponding to a = +1 and 
a = — 1 are similar. 
We get the Fig. IT51 

Estimation problems 

We must estimate the branches C\,D\, C' 1 D' 1 ,Ci,C' 1 ,C2, D 2 and 62- We know 
a priori the following : 

• The branches Cx^C'^C'^D'^ and C 2 are semi- analytic. We must check if 
they end on the abnormal direction. 
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Figure 17: 

• The branches C% , D\ , C2 and D2 are in the exp-log category and are ending 
on the abnormal direction. 

• We must compare the positions of the branches Ci , Di , C2 and D2 to 
determine which ones are in the sphere. 

All the computations are made in the general integrable case, i.e. the coef- 
ficients of the metrics a and c are analytic functions of y so that : 

a = 1 + lay + • • • 
c = 1 + 2jy H 

Our computations are based on the integral formulas (|35|l and lead to the fol- 
lowing : 
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Figure 18: a = ±1 



Proposition 4.12 (Comparison of branches C\,C2, -Da)- Let X — a^rf 
% = -fS' We ftawe i/ie estimates : 

• &ranc/i d : Z = ±A 3 + + ff (a + 7)) X 4 + o(X 4 ) 

• &ranc/i C 2 : Z = ^X 3 + opf 3 ) 

• feranc/i D 2 ; Z = ±X 3 + - f§(a + 7)) X 4 + o(X 4 ) 
and we can conclude : 

• if 7 > — a, i/ie branch C\ is in the sphere. 

• if 7 < —a, £/ie branch D2 is in the sphere. 

Remark 4.4. At the Gauss curvature of the Riemannian metric gu — adx 2 + 
cdy 2 is K = , If /3 = 0, it reduces to Hif^tll. Hence the critical value 

a + 7 = is connected to K — 0. 

If a = in the gradated form of order 0, the section reduces to : y = 0. Then 
the branch D\ does not exist (see Fig. 1171) and the branch C\ = C[ ends on the 
abnormal direction (and is in the sphere). Also the branches C' 1 D' 1 and C2 end 
on the abnormal direction, but are not in the sphere, as can easily checked. 

If a ^ 0, the branches C\, C[, C[D' 1 and C2 do not end on the abnormal 
direction. A new branch appears : D\, which is the only branch in z < that 
ends on (— r, 0) (the same is available in z > on (r, 0)). Therefore D\ is in the 
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sphere. 



Hence we know the asymptotics of the trace of the sphere with y = near 
the singularity (— r, 0) (resp. (r, 0)) in the general integrable case. Now an 
important question is to check in which class it is. In [2] it was proved that the 
sphere in the flat case is not subanalytic. Very precise evaluations of flat terms 
of branch C\ lead to the following : 

Theorem 4.13. In the general integrable case the sphere is not subanalytic. 

Remark 4.5. This result cannot be obtained by perturbation of the flat case. 
The explanation is the following. 

We proved that in the flat case the sphere is not subanalytic : 

Z = -X a - 4X 3 e~x + o(X 3 e-^) 
6 

In the general case (not only integrable) a natural idea would be to invoke some 
perturbation argument in order to check non subanalyticity. We may think that 
the previous graph is continuous with respect to the coefficients of the metrics, 
or with respect to the radius of the sphere. But this is wrong, as shown in the 
following example : 

F i = ( 1 + e v)-H- + %r-i- , *2 = ir, 9 = n 1 , 2 dx 2 + dy 2 (£<0) 



dx 2 dz dy ' (1 + ey) 



We obtain 



Z=\x" + ^X± + ■■■ + re{\ ^ 2 e 2 )lV^ + o{X^) 

This is actually not surprising, since in the step of elimination of the parameter 
k' (see |14)1. we replaced k' with its expression in function of X. But this step 
needs an exponentiation, and we know that equivalents do not pass through 
exponentiation. 

However we could expect that the expansions of X and Z in function of , k! 
(see |21 ) are continuous with respect to the coefficients. It is still wrong : 

• flat case : Z - \X 3 = -2-^ + o(^) 

. case e < : Z - \X 3 = ^An(^) - (3 + + o(fy 

Nevertheless we can observe that the analytic part of the graph is always con- 
tinuous with respect to the coefficients. Instability only appears in flat terms. 
This can be easily explained in the case e < : to compute X and Z, we need 
to evaluate some integrals. To do that, the change of variable r\ = 4-sh t is 
relevant (see an( i leads to expand X and Z as a sum of terms containing 
Argsh 2 ^f^- . Now if one wants to expand this last expression (using the formula 
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Argsh x — Ya(x + Vl + x 2 )), with k'^/X — ► 0, it is necessary to assume e fixed 
(so as r) to get : 

Argsh^^ln^ + In (l + Jl + ^A 

= ln + In2 + An(^ r 4) 

2fc'\/A ye 2 ^ / 

in order to obtain analytic expansions of X and Z , which prove that the sphere 
belongs to the log-exp category. Unfortunately in this last expression, there is 
no sense to make e — > because we needed to assume e fixed. Moreover note 
that fc' sh 2 Argsh ^t^j = + k' +o(k' ), so that this term brings new flat 
terms with coefficients having the same order as unity. 

We could now expect to have continuity with respect to parameters if we do 
not expand the Argsh's, and try to make the following reasoning : 

f . x i — > /(0, x) is not subanalytic. 

2. e i — > f(e, x) is continuous. 

Then for e ^ x i — > f(s,x) is not subanalytic. 
But this is wrong, see the following example : 

£ t 2 

/(e, t) = In t + Argsh - = In 2e H ^ + • ■ ■ : analytic in t. 

t l£ l 

/(0,£) = ln< : not subanalytic. 

So the sphere is not subanalytic. Now the main question is : in which 
category is the sphere ? In we proved that the branch C\ belongs to the 
log-exp category. A precise answer is the following : 

Proposition 4.14. We set near the singularity (— r, 0) : X = ^f, Z = 4-, 

and we have : 

• branch d : Z = An(X, X\nX, X\n 2 X, X\n 3 X, = \X 3 + ■■■ 
where An{.) is a germ at of an analytic function. Moreover the analytic 
part of Z(X) is continuous with respect to r and the coefficients of the 
metrics. 

A similar result holds for D2 ■ 

• branch D x : Z = An(y/^X, V^Xln(-A), e~^*7^r) = -^X 2 + ■■■ 
Moreover the analytic part of X(\/Z) is continuous with respect to r and 
the coefficients of the metrics. 

Corollary 4.15. In the general integrable Martinet case the sphere belongs to 
the log-exp category. 

Proof. Our estimations show that near the abnormal direction the sphere is 
log-exp. In the other directions the sphere is subanalytic, see PP. □ 
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4.2.13 Asymptotics of the sphere and wave front in the abnormal 
direction in the general gradated case of order 

We set : g = (1 + ay) 2 dx 2 + (1 + f3x + ^y) 2 dy 2 with a, f3 ^ 0. In this case the 
equation in (0, 6) obtained by projection is not integrable. In order to compute 
the asymptotics of the sphere we can use formal first integrals near the saddles. 
Moreover toric blowing-up allow us to evaluate the solution if A is fixed, see 
The technics are similar to the ones used by [SSj and others to evaluate 
the Poincare-Dulac return mapping near a polycycle for a one-parameter family 
(X £ ) of vector fields. This computation can be reduced to the evaluation of the 
Poincare-Dulac mapping near a resonant saddle : 

r(0) 

In this method we use a normal form in which : 

• the separatrices are normalized to : x' = 0, y 1 = 0. 

up to a change of time parametrization, and in some neighborhood U x £(e) C 
R 2 x R, with £(e) — ► 0, see |3H] for details. 



= ^44 a ^<o 

A 2 (e) 



= -1 



r = 



= —1 + ax{e) 



3 I 

0-2 




In our case the method has to be adapted. Indeed in our case the section 
depends on the parameter and is imposed by the geometry. In particular the 
distance of the saddle point to the section a s tends to when e — > 0. The 
method is the following. Let d be the distance to the separatrix. Then we want 
to compute : d i — ► (x(d) + r,z(d)) when d ~ (using the normal form for 
(X £ ,a e )). 

This computation generalizes the computation in the conservative case where 
d is the distance to the root of multiplicity two of the potential. 

The algorithm to evaluate step by step this application is to consider the 
fc-jet of (X e ,a E ). It is not clear a priori that the A; -jet is sufficient to compute 
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the first k terms in the expansion. However we shall prove that the 1-jet is 
sufficient to compute the first term in the expansion. It gives us the contact of 
the branch D\ with the abnormal direction. 



Proposition 4.16. Let us suppose a = (1 + ay) 2 , c = (1 + Sx + ^y) 2 with a > 

x+r 
1r ' 



Let X = , Z = -%. Then near X — the graph of the branch D\ is the 



following : 

Z = -^X 2 + o{X 2 ) 
r z a z 

Remark 4.6. Observe that in the flat case, the abnormal geodesic is not strict 
and the contact is of order 1 (see prop l4.10")l . 

Proof. The differential system is : 

dx cos 8 dy sin 8 dz y 2 cos 8 



(39) 



dt 1 + ay dt 1 + fix + jy dt 2 1 + ay 
d8 1 

— = — 7- rj- r (Xy — a cos 8 + 3 sin 8) 

dt (i + m/)(i + fix + iy y y ' ' 

Reparametrizing with : ds = y/X (± +ay )^ +j3x+iy) dt, we obtain : 

dx 1 ,„ „ . . 

— = -j={\ + Bx + ^y) cos 8 

^ = — !=(1 + ay) sin# 
ds VA 

dz 1 v 2 
_ = _ T (l + ^ + 7 y)ooBfl 

dO nr a n ft ■ „ 

— = — v Ay H = cos # = sin # 

«s VA VA 

Hence the equation governing 8 is : 

d 20 ■ „ « 2 ■ a a a P • 2 a ft ndO n 

—r + sm 8 + — sm cos 8 — sin 8 H = cos 8— = 

aV A A vA as 

Set u = 8 + tt,v = |§. Then : 
as 

dt; a 2 . ad , 2 /3 

— = sm u — sm u cos u H — — sm u H —v cos u 

ds X X s/\ 

The eigenvalues of the linearized system are solutions of : /x 2 — -^/i— (1 — ^-) = 0, 
hence: Ml = l + ^= + 0(i) , M2 = -1 + ^= + O (i). 
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Let u = Ui + V\, v = jUiiti + /i2"2- We get : 

du\ ( u\ v\ u\V\ u? u\v\ u\v\ v\ 

— = ii x u\ + O — , — , — — , -=, 



ds ^ ~ V a ' a ' a 'Va'Va'Va'Va 

dvi f u\ vf U1V1 u\ u\v\ u\v\ v\ 



^ 1 \ A ' A ' A 'VA'VA'VA'^ 
and after integration : 

, A 2 rj2 *n 4 3 

u(a) = Ae" 13 + Be? 28 + 0( — e 2 " lS , — e 2 ^ 2S , e ^ 1+ ^> -= 

V A A A VA 



e 



d3 

y^2^ e (2/ii+/i 2 )s ^^2 e (/ii+2^i 2 )s _ c 3/J2Ji 

' Va 

«(«) = ^i^ lS + H2Be» 2S + 0(- • • ) 
where ^4 and S are constants to determine. 

The section is y = 0, hence : v = cosu + sinu = + -^=u + O (^j) 

Let be the parameter corresponding to the hnal time t — r, i.e. : y(0) = 
y(sf) = 0. Putting these conditions in the previous equations we obtain : 

Hence : 

9(a) + 7T = u(a) = I^e" 1 ('-"') + -^=e^ s + O (^e/*i'-'/\ \e" 2S ) (40) 
vA vA \ A A / 

To get y, just note that : y = — + f cos 6* — sin#, hence : 

y(s) = "X + f eMl(S " S/) + f e " 2S + (3^72) (41) 

Then we have to compute x, which amounts to integrating equation (|39|) . 
We get : 

l + M S )=e-^ + ^-fe-^ + 0(J 7I ) (42) 

The computation of z is then similar and we obtain : 

a 2 l-e"^ s / 1 \ 
<») = -TXi H +° T2 ( 43 ) 



2A 2 /3 VA 2 

It remains to estimate Sf. From the equation : = ^(1 + a V) (1 + fl x + 72/) 
we get : 

This leads to the conclusion : Z = -p^jX 2 + o{X 2 ) 

□ 
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Remark 4.7. Another way to compute this expansion is to use the theory devel- 
oped in 05| j which states that the so-called I/°°-sector has a contact of order 2 
with the abnormal direction, and moreover gives an explicit formula to estimate 
the contact. 

The previous method cannot be applied to study the contact of branches C\ 
and D2 with the abnormal direction, because in the phase plane of the pendu- 
lum these branches correspond to a global computation of return mapping, and 
thus the calculations cannot be localized near a saddle as previously. Anyway 
inspecting carefully the system leads to the following : 

Lemma 4.17. In the general gradated case of order the contact of branches 
C\ and D2 with the abnormal direction is : 

Z={\ + 0(r))X 3 + o(X 3 ) 

Note that contacts are still in the polynomial category. 

Proof. We have : y = v/y/c, hence ||y||oo = 0(r). On the other part : x = 

2 

u/ (1 + ay), and thus : x = u(l + O(r)). In the same way : z — + O(r)). 

Then the result in the flat case leads easily to the conclusion. □ 

Remark 4.8. From our previous study we can assert that minimizing controls 
steering to points of C\ (resp. D2) are close to the abnormal reference control 
in L 2 -topology, but not in I/°°-topology. It is a crucial difference with the branch 
Dl 

Concerning the transcendance of this branch D\, the following fact was 
proved in 02j : 

Proposition 4.18. In the general gradated case of order 0, if a 7^ then the 
branch D\ is C°° and is not subanalytic at x = — r, z = 0. 

Corollary 4.19. In the general gradated case of order 0, if the abnormal min- 
imizer is strict then the SR spheres with small radii are not subanalytic. 

Proof. Let A — (— r, 0, 0) denote the end-point of the abnormal trajectory. We 
shall prove that D\ is not subanalytic at A. The method is the following. First 
of all the Maximum Principle gives a parametrization of minimizing trajectories 
steering to points of D\ . Then we prove that the set of Lagrange multipliers 
associated to these points (i.e. end-points of the corresponding adjoint vectors) 
is not subanalytic. Finally we conclude using the fact that, roughly speaking, 
these vectors coincide with the gradient of the sub-Riemannian distance (where 
it is well-defined). These facts are summarized in the following : 

Lemma 4.20. To each point q of D\ is associated a control u, and we denote 
by (ip(q),ip°(q)) an associated Lagrange multiplier. Then we set : 



/96O1 
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where ip x (resp. ip z ) is the projection on the axis x (resp. on the axis z) of the 
vector ip. If the set C is not subanalytic then the curve D\ is not subanalytic. 

Proof of the Lemma. Let (9(t))o^ t <i be a parametrization of the curve D\ such 
that q(0) = A. For each r let u T be a control such that E(u T ) — q(r), and let 
(V'tj'0?) be an associated Lagrange multiplier, i.e. : 

i) T AE{u r ) = -ip°dC(u T ) 

Then : ^ T .£g(r) = ^j T .dE(u T ).£u T = -^° T dC(u T ).£u T = -f^CK). 
Moreover for each r the point q(r) belongs to the sphere S(0, r), hence C(u T ) = 
r, and thus : tj) T .4-q{T) = 0. Therefore in the plane (y = 0) the vectors of 
the set £ are unitary normal vectors to the curve D 1 . Then the conclusion is 
immediate. □ 

With notations of Proof of Proposition 14.161 we are now lead to study a 
family of vector fields (X £ ) depending on the parameter e = -7=, in the neigh- 
borhood of a saddle point u = v — 0. For the section E corresponding to y = 
we estimate the return time, i.e. the time needed to a trajectory starting from E 
to reach again E ; then we claim that this time is t = r. This gives us a relation 
between 8{r) and A, thus between p x (r) and p z {r). Then one has to show that 
this relation is not subanalytic. We proceed in the following way. First of all 
recall that 

l + (3x{s) = e ~^ s + 0(j), y(s)=0(j) 

We need a result which is independant of the parameter e = . So it is no use 
trying to write an analytic normal form, since the saddle may be resonant. On 
the other hand C k normal forms (see are not enough because flat terms that 
we aim to exhibit disappear up to a 0(e k ). However near the saddle separatrices 
of X £ are analytic in u, v, e, and actually there exists an analytic change of 
coordinates (ui,Vi) = An(u,v) (in the sequel An(.) denotes an analytic germ 
at 0) such that in these new coordinates separatrices are u\ — 0, v± — 0, and 
the system is : 

where ^ 2 (-^) are the eigenvalues of the saddle ; in particular : fj, 2 (^=) = 

— 1 + -A= + 0(I). Moreover we have : u — u± +«i +o(^=), v = u\ — v\ +o(-^=), 
therefore the section is E : v\ = u\ + + o(-^=). Let Sf denote the pa- 
rameter corresponding to the return time, i.e. (ui(s f) , v±(s t)) G E. We have : 

»'<°> = 7T + °<7x)- T1 - : 

Jo Jvt(0) A*2(^K(l + o(^)) ^2 VA M0) 



53 



On the other part : £ = -±(l + ay)(l + /3x+~/y) = -^e 7x8 + 0(^). Hence 
we get : 

f' f dt. l-e~vr Sf 1. 
And thus : s f = -VA ln (1 ~^ r) +0(4-). Putting into (gHJl we obtain finally : 

In particular «i (s/) is not an analytic function in . 

We know that Vi(sf) — An(u(sf),v(sf)) ~ u ( s /)~"( 8 /) _ Moreover, on the 
section E, we have : v(sf) = —^cosu(sf) + -^=sinu(s/). Hence : vi(sf) = 
An(u(sf), = g + • • ■ • From the Implicit Function Theorem in the ana- 
lytic class we get : u(sf) — An(vx(sf), Therefore u(sj) is not an analytic 
function in -^=, So the set C is not subanalytic, which ends the proof. □ 

Unfortunately we have no general result similar to Corollary 14. 151 We think 
that the log-exp category is not wide enough in the non integrable case. Indeed 
due to the dissipation phenomenon observed in the pendulum representation if 
7^ we cannot expect to keep the analytic properties required in the definition 
of log-exp functions. Moreover in the phase plane of the pendulum, the foliation 
is not a priori integrable in the analytic category for any value of the parameter. 
We can observe that if we fix A to 1 and evaluate the Poincare-Dulac mapping, 
it is pfaffian if and only if X e (A = 1) is C^-integrable (see [34"1 135] ) . 

Hence we conjecture : 

Conjecture 4.21. If (3 ^ then SR Martinet spheres are not log-exp, even not 
pfaffian. 

Hence we should try to extend the log-exp category to a wider category in 
which analyticity would be replaced with some asymptotic properties. That is 
why we should be interested in Il'Yashenko's class of functions (see \22\). In- 
deed several years ago the Dulac's problem of finiteness of limit cycles was solved 
independantly by Ecalle and Il'Yashenko ; in his proof, Il'Yashenko introduces 
a very wide class of non-oscillating functions to describe Poincare return map- 
pings. Actually he needs to expand real germs of functions into terms having 
not only the order of x n but also into flat terms. His category is contructed by 
recurrence as follows. Let Mq be a class of functions that can be expanded in 
an unique way into an ordinary Dulac's series, i.e. series of the form : 



So = cx Uo +^2 Pi (In x)x^ 



where c > 0, the Pi's are polynomial and {vi) is an increasing sequence of posi- 
tive numbers going to infinity. We only sketch the first step of the recurrence. 
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By definition germs of the class Mi can be expanded in series containing flat ex- 
ponential terms whose coefficients belong to Mq. For instance a super- accurate 
series of some germ / may be of the form : 

oo 

Si = a (x) + 2J ai(x)e~^ 

i=l 

where € Mq and is an increasing sequence of positive numbers tending to 
infinity. It is a generalization of ordinary series in so much as the usual Dulac's 
series of / is a (x) . 

The complete definition of super-accurate series then goes by recurrence (see 
\22\). Their interest is all in the fact that the application / > f is one-to-one, 
where / is the super-accurate series associated to /. 

The relation with our problem is the following. We deal actually with 
Poincare return mappings in the phase plane of the pendulum, and their study 
is crucial to estimate the spheres. The difference is that our pendulum depends 
on parameters (namely A) ; hence our problem is related to Dulac's problem 
with parameters, i.e. the Hilbert's 16th problem. 

Hence we should try to construct a category of functions similar to the one 
introduced by IPYashcnko, but with parameters. In any case it is a possible 
way to try to solve the problem of transcendence in SR geometry. 

Conjecture 4.22. In the general non integrable case SR spheres belong to 
some extended Il'Yashenko's category. 



4.2.14 Conjecture about the cut-locus : the Martinet sphere in the 
Liu-Sussmann example 

We shall construct the cut-locus in the Liu-Sussmann example |29j . the reason- 
ing being generalizable to compute the generic SR-Martinet sphere. The model 
is the following : 

D = Ker u , uj = (1 + ey)dz - —dx , g = — ^ + dy 2 

The model is non generic because it is conservative ; moreover the Lie algebra 
generated by the orthonormal frame is nilpotent. In the cylindric coordinates 
the geodesies equations are : 

V 2 

x = (1 + sy) cos 8 , y = sin 6 , z = — cos 

= -{Px£+PzV) , P3 = A 

and the angle evolution is the pendulum : 6 + Asin# = if A 7^ 0. Using 
the symmetry : (x,y,z) 1 — ► (— x, y, — z) we can assume A > 0. The abnormal 
geodesic starting from : t 1 — > (±f, 0, 0), is strict if and only if e 7^ 0. We 
may assume s ^ 0. Introducing s = i\/A, and denoting by ' the derivative with 
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respect to s, the pendulum is normalized to : 9 + sin 8 = 0. The constraint 
y = defines the section £ : 8 = —p x e, which can be written : 8' = — EC ^L 9 . 

The geodesies corresponding to A = are globally optimal if the length r 
is small enough. They divide the sphere S(0, r) into two hemispheres and we 
compute the cut-locus in the northern hemisphere (A > 0). 




Figure 19: 

If e = 0, the SR-sphere is the Martinet flat sphere. The abnormal line is not 
strict and cuts the equator A = in two points. The cut-locus is the plane y = 
minus the abnormal line, in which, due to the symmetry (x, y, z) i — ► [x, —y, z), 
two normal geodesies are intersecting with the same length. It is represented 
on Fig. HI 

When e ^ 0, the section and the pendulum are represented on Fig. [501 




Figure 20: 

(7 is given by 8' = — £C ^ e ■ When A — > +oo, the section tends to 8' — 0, and 
when A — ► the points M and m tend to oo. 
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We can easily construct the cut-locus on the small sphere using the following 
conjectures : 

• Only the geodesies where the section a is in the configuration of Fig. [2U] 
have cut-points (this is justified by the fact that when A = the geodesies 
are globally optimal if the length is small enough). 

• A separatrix has no cut-point. 

The trace of the sphere with the Martinet plane has been computed in ^U] . 
An important property is that the rotating trajectories of the pendulum near 
the separatrix have a cut-point located in the plane y — 0, corresponding to its 
second intersection with y = 0. We represent on Fig. [2]the construction of the 
cut-locus. 



The cut-locus has two branches Lc and Ljj corresponding respectively to 
oscillating and rotating trajectories. They ramify on the abnormal direction 
A which is not contained in the cut-locus. The extremities of the branches 
Lc and Ljj are conjugate points corresponding respectively to 8(0) = it and 
8(0) = 0. The branch Lc has only one intersection with y = which corresponds 
approximatively to 8(0) = \. The cut-locus is not subanalytic at A but belongs 
to the log-exp category, and thus from |2()j : 

Lemma 4.23. The cut-locus is C°° -stratifiable. 

To generalize this analysis we must observe the following. The respective 
positions of the branches C\ , Di are given in Section 14.2.121 Here the curve 
Z?2 is above and hence the rotating trajectories are optimal up to the second 
intersection. Also the integrability of the geodesic flow is not crucial and in 
general the branch is not contained in the plane y — 0. We make the 
following conjecture. 

Conjecture 4.24. 1. In the Martinet case the cut-locus is C 1 -stratifiablc. 




A = 



L D 



Figure 21: 
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2. In the generic case the cut-locus has two branches in the northern hemi- 
sphere ramifying at the end-point of the abnormal direction. 

5 Some extensions of Martinet SR geometry and 
microlocal analysis of the singularity of the SR 
sphere in the abnormal direction 

5.1 Non properness and Tangency Theorem 

This analysis is based on the sub-Ricmannian Martinet case, where it was shown 
in the previous Section that the exponential mapping is not proper and that in 
the generic case the sphere is tangent to the abnormal direction. This fact is 
actually general and we have the following results (see |40j ) . 

Consider a smooth sub-Riemannian structure (M, A, g) where M is a Rie- 
mannian n-dimensional manifold, n 3, A is a rank m distribution on M, and 
g is a metric on A. Let go G M ; our point of view is local and we can assume 
that M = R" and go = 0. Suppose there exists a strict (in the sense of definition 
I3.2fl abnormal trajectory 7 passing through 0. Up to reparametrizing we can 
assume that A = Span {F±, . . . , F m } where the system of F^s is <?-orthonormal. 
Then the sub-Riemannian problem is equivalent to the time-optimal problem 
for the system : 

m 

g = £«i-Pi(g), q(0) = (46) 

i=l 

where the controls satisfy the constraint X)I™=i u i ^ Suppose further that 7 
is associated to an unique strictly abnormal control. Then : 

Theorem 5.1. The exponential mapping is not proper near 7. 

Proof. Set A = j(r), r > 0. Let (A n ) be a sequence of end-points of minimizing 
normal geodesies q n converging to A. To each geodesic q n is associated a control 
u n and an adjoint vector (p n ,Pn)- As q n is normal we may suppose that p° = 
— 5. Let ip n the end-point of the adjoint vector p n . Then if E denotes the 
end-point mapping and C is the cost (here the cost is quadratic in the control), 
we have the following Lagrange multiplier equality : 

^ n .dE(u n ) = ^dC(u n ) 

If the sequence ip n were bounded then up to a subsequence it would converge to 
some ip £ M n . Now since the u n are minimizing the sequence (u n ) is bounded 
in L 2 , hence up to a subsequence it converges weakly to some u £ L 2 . Using the 
regularity properties of the end-point mapping (see |4Up. we can pass through 
the limit in the previous equality and we get : 

i).dE(u) = ^dC{u) 
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and on the other part : A = E(u). It is not difficult to see that u has to be 
minimizing, and then we get a contradiction with the fact that 7 is strict. □ 

Remark 5.1. Conversely if the exponential mapping is not proper then actually 
there exists an abnormal minimizer. This shows the interaction between abnor- 
mal and normal minimizers. In a sense normal extremals recognize abnormal 
extremals. This phenomenon of non-properness is characteristic for abnormal- 
ity. 

This non-properness is actually responsible for a phenomenon of tangency 
described in the following Theorem (see |42| for a more general statement): 

Theorem 5.2. Consider the SR system \4o\j and suppose there exists a min- 
imizing geodesic 7 associated to an unique strictly abnormal control u. Let 
A G S(0, r) be the end-point 0/7. Assume (c(t))o<t^i * s a C 1 curve on 5(0, r) 
such that lim <t(t) = A. Then : lira a (r) € Im dE(u). 

In particular if S(0,r) is C 1 -stratifiable near A then the strata of S(Q,r) 
are tangent at A to the hyperplane Im dE(u) (see Fig. 12211 . Moreover if B 
is a C 1 -branch of the cut-locus ramifying at A then B is tangent at A to this 
hyperplane. 




Figure 22: Tangency phenomenon 



Proof. For each r the point <t(t) is the end-point of a minimizing geodesic, 
and we denote by (p T ,p°) (resp. u T ) an associated adjoint vector (resp. an 
associated control). Let (ip r , ip®) be the end-point of this adjoint vector. We 
may suppose that it is unitary in R" x M. We have : 

TP T .dE{u T ) = -ij)°dC{u T ) 
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Using the same reasoning as in the Proof of Theorem l5.1l we get that \\ip®\ \ ~ * 
as t — > 0. To conclude it suffices to show that ip T is normal to the curve <t(t). 
Indeed the previous equality implies : 

dr dr 
But C(u T ) is constant (equal to r) and E(u T ) — a(u T ), and thus : 

dr 

which ends the proof. □ 



5.2 The tangential case 
5.2.1 Preliminaries 

In this Subsection we shall make a brief analysis of the so-called tangential case. 
According to Section 2.2.2 the distribution D = Ker to can be reduced to 
one of the normal forms : 

• elliptic case : uj e = dy — (exy + ^- + xz 2 + mx 3 z 2 )dz 

• hyperbolic case : u>t = dy — (exy + x 2 z + mx 3 z 2 )dz 

where s = ±1. The parameter e is a deformation parameter whose intro- 
duction will be justified later. 

A general metric g is then defined by : a(q)dx 2 + 2b(q)dxdz + c(q)dz 2 where 
a, 6, c can be taken as constant in the nilpotent approximation of order —1. Our 
study is far to be complete and we shall describe briefly the case g = dx 2 + dz 2 . 

The general case of order —1 depends on a parameter A and corresponds to 
a 6-dimensional nilpotent Lie algebra. It contains both elliptic and hyperbolic 
cases. It is the Lie algebra generated by Fi, F2 with the following Lie brackets 
relations : 

Fs = [Fi,F2\ , F4 — [F3,Fi] 
F5 = [F3,F 2 ] , F 6 — [F 4 ,Fi] 
[F5,F 2 ] — XF 6 , 

and all other Lie brackets are 0. 

Introducing Pj =< p,Fi > the geodesic equations are given by : 

A = P3P2 , Pi = -P3P1 , A = P4P1 + P5P2 , 

P 4 = P 6 Pi , A = AP 6 P 2 

and Pq is a Casimir first integral. The value A = represents the bifurcation 
between the two cases. 
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5.2.2 Abnormal geodesies 



Elliptic case The abnormal geodesies are contained in the Martinet surface : 

ey + x 2 + z 2 + 3mx 2 z 2 = 
and are solutions of the equations : 

x = (2z + 6mx 2 z) — e(— - — h 2mx 3 z 2 ) 
3 

i = — (2x + 6mx 2 z) 

From |45| . the singularity x = z = is a weak focus and a spiral passing 
through is with infinite length. Since any minimizer is smooth no piece of 
abnormal geodesic is a minimizer when computing the distance to 0. Using the 
general result of [Q, the sphere of small radius is the image by the exponential 
mapping of a compact set and it is subanalytic. This is also clearly shown by 
numerical simulations and the sphere is represented on Fig. 1231 



Figure 23: Elliptic case 

By taking e — 0, the Martinet surface becomes : x 2 + z 2 + imx 2 z 2 = 
and reduces near to : x = z = 0. Hence the spiral disappears. Since the 
weight of x, z is one and the weight of y is four, it corresponds to the nilpotent 
approximation of order —1 where m is 0. 

Hyperbolic case The Martinet surface is given by the equation : 

ey + 2xz + 3mx 2 z 2 = 
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and the abnormal geodesies are solutions of : 

x = 2x — x 2 z(e — 6m) — 2mx 3 z 2 
z = — (2z + 6mx 2 z 2 ) 

and the singularity at x = z = is a saddle point. The two lines x — and z = 
are optimal for the metric dx 2 + dz 2 . Hence they play a role when computing 
the distance to 0. Numerical simulations show that the sphere is not the image 
of a compact set. This can be seen on Fig. because the sphere cannot be 
numerically represented in the abnormal direction (there is a hole). It is similar 
to the situation encountered in the Martinet case. The sphere is pinched in both 
abnormal directions. 

The nilpotent approximation of order —1 is obtained by taking e = and 
r7i = 0. The Martinet surface becomes : xz = and the two lines x = and 
z = remain abnormal geodesies. 




Figure 24: Hyperbolic case 



5.2.3 Normal geodesies 
Elliptic case We take the frame : 

Fl = lk ' F2 = lfz + ( £XV +Y +Xz2 + mx3z2 )^~ ' 
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and we introduce : p =< p, Fi(q) >. The geodesies equations are : 
* = Pi 

y = P2(exy + — + xz 2 + mx 3 z 2 ) 
z = P 2 

Pi = -(ey + x 2 + z 2 + 3mx 2 z 2 )P 2 P 3 
P2 = (ey + x 2 + z 2 + 3mx 2 z 2 )P 1 P 3 
P 3 = -exP 3 

and they can be truncated at order —1 by making e = m = 0. In this case P 3 is 
a first integral and we can set : P 3 = A. Moreover if we introduce : Pi = sin#, 
P2 = cos 8, the equations become : 

x = sin 9 , i = cos 6 , y = cos 6 ( — + xz 2 ) 

9 = -(x 2 + z 2 )P 3 , P 3 = X 
They can be projected onto the space (x, y, 9) and the foliation is defined by : 

x = sin 9 , i = cos 9 , 9 = —(x 2 + z 2 )\ 

It is not Liouvillc-integrable but the equations can be integrated by quadra- 
tures, see |36| . Using polar coordinates : x — r cos ip,z = rsin-0, it becomes : 

• cos(6> + VQ . • (a 1 M a \ 2 

x = , r — smw + w) , 9 = — Ar 

r 

The important property is the following : 

Lemma 5.3. The sign of 9 is constant and reparametrizing the equation can be 
rewritten : 9 — 1 as in the contact case. 

Numerical simulations 

The geodesies equations can be integrated numerically. The projections in 
the plane (x, z) of the geodesies starting from are flowers with three petals 
(they are circles in the Heisenberg case), see Fig. |53J 
Hyperbolic case We take the frame : 

Fi = ^-, F 2 = ^- + {exy + x 2 z + mx 3 z 2 )-^- , F 3 = ^~ 
ox oz ay ay 

and we introduce : Pi —< p,Fi(q) >. The geodesies equations are : 
x = Pi, y = P2{exy + x 2 z + mx 3 z 2 ), z = P 2 
Pi = -(ey + 2xz + 3mx 2 z 2 )P 2 P 3 
P 2 = {ey + 2xz + 3mx 2 z 2 )PiP 3 
P 3 = -exP 2 P 3 
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Figure 25: 

and they can be truncated at order —1 by making e = m = 0. In this case P3 is 
a first integral and we can set : P3 — A. Moreover if we introduce : P\ = sin9, 
P2 = cos 8, the equations become : 



x = sin f 



y — cos a x z 



6 = -2xzP 3 , P 3 = A 
They can be projected onto the space (x, z, 6). 





z z 






X 


X 



Figure 26: 



Numerical simulations 

The projected equations : x = smO, z = cos 9,9 = —2xzX can be integrated 
numerically and the solutions compared with the pendulum : x = sin 9,i — 
cos 9, 9 — —Ax, see Fig. EH1 The behaviour is quite chaotic and 9 exhibits oscil- 
lating and dissipative phenomena. Also it shows in the plane (x, z) a coupling 
effect between the two abnormal directions. 

The SR sphere in the tangential hyperbolic sphere is represented on Fig. [23 
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Figure 27: Tangential hyperbolic sphere 



5.2.4 Conclusion 

The elliptic situation is similar to the contact case. The analysis of the hy- 
perbolic case is intricate. A tool to study the sphere is to introduce as in the 
Martinet case a return mapping by taking the intersections of the geodesies with 
one of the planes : x — or z — 0. The non properness of this application can 
be checked numerically. 



5.3 The Engel case and left-invariant SR geometry on 
nilpotent Lie groups 

If q = (x, y, z, w), we consider the system in R 4 : 



d d y 2 d d 
1 dx ^ dz 2 dw' 2 dy 



We have the following relations : F3 = [Fi, F2] 



[[F 1 ,F 2 ],F 2 ] 



■J^ and [[Fi, F 2 ], Fi] — 0. Moreover all Lie brackets with length greater than 4 
are equal to zero. Set : 



Li = 



( 





\ 





1 



Lo = 







0/ 



\ 





\ J 



and define the following representation : p(F\) = Li,p(F 2 ) = L 2 which allows 
to identify the previous system in 1R 4 to the left-invariant system R = (uiLi + 
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u 2 L 2 )R on the Engel group G e , here represented by the nilpotent matrices : 

/I 92 93 94 \ 

o i 9i i 

1 qi 
V 1 / 

The weight of x, y is one, the weight of z is two, and the weight of w is three. 
For any sub-Riemannian metric on G ei the approximation of order —1 is the 
flat metric g = dx 2 + dy 2 . Any sub-Riemannian Martinet metric can be written 
g = adx 2 + cdy 2 and can be lifted on G e . 



5.3.1 Paramatrization of geodesies in the flat case 

Non trivial abnormal extremals are solutions of : 

Pl = P2 = {Pi,P 2 } = 0, u 1 {{P 1 ,P 2 },P 1 } + u 2 {{Pi,P 2 },P 2 } = 
Set p = (p x ,p y ,p z ,p w ). We get : 

y 2 

Px + PzV + Pwy = Py = P*> + VP™ = PwU 2 = o 
This implies p w ^ and thus U2 = 0. The abnormal flow is given by : 

y 2 

x = ux, y = 0, z = uxy, w = m — 

where |mi| = 1 if the parameter is the length. 

To compute normal extremals, we set Pi =< p,Fi(q) >,i = 1,2,3,4 and 
H n = \{P 2 + P 2 ) ; we get: 

Pi=P 2 P 3 , P 2 = -PiP s , P^=P 2 Pi, P 4 = o 

Parametrizing by the length H n = i, we may set : Pi — cos9, P 2 = sin#, and if 
^ kir we get : 8 — —P3, = —P 2 P±. Denote P4 — A, then this is equivalent 
to the pendulum equation : 

6 + X sin 6 = 

Let L denote the abnormal line starting from : £ 1 — > (±i, 0, 0, 0). It is not 
strict and projects onto 9 = kir. 

In order to obtain an uniform representation of normal geodesies, we shall use 
the Weierstrass elliptic function V . Indeed the system admits three integrals : 
P\ + p 2 = 1) and two Casimir functions : -2PiP 4 + P 3 2 = C ct P 4 = A. Using 
Pi = P 2 P 3 we get : A = -CP\ - i\P 2 + A, which is equivalent, with A^0 

and A ^ 0, to the equation : Pi = -2A(F 1 3 + %P\ -Pi+D). Let V[u) denote 
the Weierstrass elliptic function (cf solution of : 

V'{u) = -2^{V{u) - ei)(T(u) - e 2 )(T(u) - e 3 ) 
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where the complex numbers Ci satisfy e\ + e 2 + e$ = 0. Set g% = — 4(e2e3 + 
e3ei + e\e?) and g$ = 4eie2e3 ; then : V(u) = 4"P 3 (w) — giPiu) — g 3 . The 
function 'P(u) can be expanded at in the following way : 

v{u) = h + h q2u2 + h q " uA + --- 

Hence the solution Pi can be written : aP{u) + b. Then we can compute P 2 and 
Pz using the integrals, and x, y, z, w can be computed by quadratures. We find 
again oscillating and rotating solutions of the pendulum using Jacobi elliptic 
functions given by the formulas : 

(V{u)-e x \^ (V(u)-e 2 \ i 

cn!i= —7 , , dn u = —. , 

\V{u)-e 2 ) \V(u)-e 3 J 

5.3.2 Heisenberg and Martinet flat cases deduced from the Engel 
case. Blowing-up in lines 

Note that the two vector fields and commute with F± and Fo ■ The Engel 

oz aw 

case contains the flat contact case and the flat Martinet case which are given 
by the following operations : 

• Setting p z = 0, we obtain the geodesies of the Heisenberg case. 

• Setting p w — 0, we obtain the geodesies of the Martinet flat case. 

The interpretation is the following. 

Lemma 5.4. We obtain the Martinet flat case (resp. Heisenberg) by minimizing 
the SR distance to the line {Oz) (resp. (Ow)). 

Indeed the condition p z = (resp. p w = 0) corresponds to the transversality 
condition. It may be observed that, since the SR distance to a line is more 
regular than to a point, the SR distance in the Engel case has at least all 
singularities of the Heisenberg and Martinet flat cases. 

Another way to get the Martinet flat case is to use the following general fact 
from [7| : 

Lemma 5.5. The Martinet flat case is isometric to {G e /n, dx 1 + dy 2 ) where H 
is the following sub-group of G e : {exp t[Fi,F2] / t £ M}. 

The Engel case can be imbedded in any dimension, for more details about 
left-invariant SR-geometry on nilpotent Lie groups see |39j. 

In the Martinet flat case, using the uniform parametrization of the geodesies 
by elliptic functions, the sphere is evaluated in any direction. 

Next we give a description of the SR sphere in an abnormal direction when 
the flag associated to the distribution D satisfies D 3 ^ D 2 , see Section 3.4. 
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5.4 Microlocal analysis of the singularity of the SR sphere 
in the abnormal direction 

The aim of this Section is to stratify the singularity of the SR sphere in the ab- 
normal direction. We use symplectic geometry. This leads to a stratification of 
the solutions of the Hamilton-Jacobi-Bellman equation viewed in the cotangent 
bundle. 

5.4.1 Lagrangian manifolds and generating mapping 

Definition 5.1. Let (M,oj) be a smooth symplectic manifold and L C M be a 
smooth regular submanifold. We say that L is isotropic if the restriction of u> 
to TL is equal to zero, and if dim L = ^dim M then L is called Lagrangian. 

The following result is crucial, see |32| . 

Proposition 5.6. Let {M,uj) be a In- dimensional manifold and L C M be 
a Lagrangian submanifold. Then there exists Darboux local coordinates (q,p) 
and a smooth function S(qi,pj) where I = {1, . . . ,m\, I — {m + 1, . . . ,n} is a 
partition of {1, . . . , n} such that L is given locally by the equations : 

dS dS 

P* = ' q i = 

oqi opj 

Definition 5.2. The mapping S which represents locally L is called the gener- 
ating mapping of L. 

Definition 5.3. Let L be a Lagrangian manifold and ff the standard projec- 
tion (q,p) i— > q from TM onto M. The caustic is the projection on M of the 
singularities of (£,IT). 

5.4.2 Lagrangian manifolds and SR normal case 

Consider the SR problem : 

q = U\F x (q) + u 2 F 2 (q), q S M 
where the length of q is : 

L(q)= [ (ul(t) + u 2 2 (t))idt 
Jo 

We use the notations of Section 3.2. Let Pi =< p, Fi(q) >, i = 1, 2, the Hamil- 
tonian associated to normal geodesies is given by H n = \{Pf + P 2 2 )- 

Let 1 1— > j(t),t £ [0, T] be a reference one-to-one normal geodesic. We assume 
the following : 

Hypothesis : We assume that the reference geodesic is strict, 
i.e. there exists an unique lifting [7] of 7 in the projective bundle 
P{T*M). 



08 



Notations 

• exp 7 ( ) is the exponential mapping. If the geodesies are parametrized by 
arc-length H n — i, it is defined by t ^ U(q(tj) where t *— ► q(t) is a solution 
of H n starting from 7(0) at time t = 0. 

• L t = cxp tH n (T*^M), where exp tH n is the local one-parameter group 
associated to H n . 

The length of a curve does not depend on the parametrization and the op- 
timal control problem is parametric. This induces a symmetry which has to be 
taken into account when writing Hamilton- Jacobi equation in the normal case. 
Indeed we have : 

Lemma 5.7. The solutions of H n satisfy the relation : 

q(t,q 1 ,Xp 1 ) = q(\t, qi ,pi), p(t,q 1 ,Xp 1 ) = Xp(Xt,q 1 ,p 1 ) 
The following results are standard : 

Proposition 5.8. 1. L a = T*^M is a linear Lagrangian manifold, and for 
each t > 0, L t is a Lagrangian manifold. 

2. The time t c is conjugate along 7 if and only if the projection H : L tc — > M 
is singular at 7(t c ). 

3. Assume that geodesies are parametrized by arc-length t, and let 

W= |J exptH n (T; i0) Mn(H n = ±)) 

0<t^T 

where T < t\ c (first conjugate time along 7. Then £ = Tt{W) is a central 
field along 7. 

Remark 5.2. • The caustic of L t is the set of conjugate points which can be 
analyzed using Lagrangian singularities. 

• We represent locally W by an Hamilton- Jacobi or wave function defined as 
follows. We integrate the normal flow starting from 7(0) and parametrized 
by arc-length : Pi + P| = 1. By setting Pl(0) = cos 6, this gives us the 
family of geodesies : 

E :p= (0,Ai,...,A„_2,t) e §1 x JR"" 1 ^ M 

and p is eliminated by solving the equation E{p) = q near j(t) using the 
Implicit Function Theorem. Beyond the computations need the Prepara- 
tion Theorem and Legendrian singularity theory. 

Next we describe the tangent space to the Lagrangian manifold. 
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Definition 5.4. We denote by (V n ) the variational equation : 

F) I i 

8§=?pW))8Q (47) 

along the reference trajectory t i— » j(t). This Hamiltonian linear equation is 
called Jacobi equation. A Jacobi field J(t) = (6q(t), 5p(t)) is a nontrivial solution 
of gZJ- It is called vertical if 8q(0) = 0. 

Proposition 5.9. 1. Let L t — exp tH n (T*^M). Then the space of vertical 
Jacobi fields is the tangent space to L t for t > 0. 

2. Assume we are in the analytic category. Let J(.) be a vertical Jacobi field 
and let e i— > a(e) be an analytic curve such that a(0) — J(0). If Y is an 
analytic vector field on T*M such that Y(f)(fS)) = a(0), then t i— * J(t) is 
given for t small by the Baker-Campbell-Hausdorff formula : 

J(t) = ^2~ i ad n H n (Y)m)) 
^-^ n! 

A consequence of Lemma 15.71 is : 

Lemma 5.10. Let 7(0) = (qo,po) and consider the curve a(e) = (qo,pa + epo). 
Then it is a vertical curve, and if Ji is the associated Jacobi field then Tl(Ji(t)) = 
tj(t). 

5.4.3 Isotropic manifolds and SR abnormal case 

Consider the system q = u\Fi(q) + ^2^2(9)- According to Section 3, the abnor- 
mal geodesies are solutions of the equations : 

dq _ dH a dp _ dH a 
dt dp ' dt dq 

where H a = U\P\ + ^2^2- They are contained in : 

P 1 =P 2 ={P 1 ,P 2 } = 

and the abnormal controls are computed using : 

ui{{Pi, P 2 }, Pi} + u 2 {{P u P 2 }, P2} = 

Assumptions Let t 1— > j{t),t e [— T, T] be a one-to-one abnormal reference 
geodesic. One may assume that it corresponds to the control u 2 — 0. We 
suppose that the following conditions are satisfied along 7 for the couple F 2 ), 
see Section 2.5. 



(Hi) The hrst order Pontryagin's cone K(t) = Span {ad k F\.F 2 \ 1 /fee N} has 

\\~t(ty 



codimension one and is generated by {F 2 , . . . , ad n 2 Fi.F 2 }\ 



70 



(H 2 ) If n ^ 3, for each t, Fi(j(t)) <£ Span {ad k F 1 .F 2h / k = . . . n - 3}. 

(#s) {P2,{ J Pi,P 2 }}^0along 7 . 

Notations 

• Under the previous assumptions 7 admits an unique lifting [7] = (7,p 7 ) in 
P(T*M). One may identify locally M to a neighborhood J7 of 7 (0) = in 
IR™. Let U be a neighborhood of p 1 in P(TqU). We can choose V small 
enough such that all abnormal geodesies starting from {0} x V satisfy 
the assumptions (Hi — H3). We denote by S r the sector of U covered 
by abnormal geodesies with length ^ r and starting from {0} x V. This 
defines a mapping denoted Exp. The construction is represented on Fig. 




P(T*U) 



Figure 28: 



• On n = T*M\({{P 1 ,P 2 },P 2 } = 0) let H a be the Hamiltonian H a = 
Pi + uP 2 where u = — {{p^p^'p^} , and let exp tH a be the one-parameter 
local group. We denote by LI (resp. I t a ) the image of T*^M (resp. 

T; {0) Mn(p 1 = P 2 = {p 1 ,p 2 } = o)). 

Lemma 5.11. On Q, L® is a Lagrangian submanifold and I® is an isotropic 
submanifold. 

5.4.4 The smooth abnormal sector of the SR sphere 

Lemma 5.12. Consider the SR problem where 7 is an abnormal reference tra- 
jectory and assume (Hi — H%). It can be identified to a trajectory of Fi where 
the system (Fi,F 2 ) is orthonormal. Then the abnormal geodesic is strict, and 
there exists r > such that if the length of 7 is less than r then 7 is a global 
minimizer. 

Proof. Under assumption (Hi) the first order Pontryagin's cone along 7 has 
codimension one, and from (H3) 7 is not a normal geodesic. The optimality 
assertion follows from jS], see also - d 

Hence the end-point of 7 belongs to the sphere. Moreover r can be estimated 
and the estimate is uniform for each abnormal geodesic C 1 -close to 7. Therefore 
we have : 
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Proposition 5.13. For r small enough S r is a sector of the SR ball homeo- 
morphic to CU-C, where C is a positive cone of dimension n — 3 if n 4 and 
1 ifn = 3. Its intersection with the sphere consists of two (C°° or C a ') surfaces 
of dimension n — 4 if n 4 and reduced to two points if n = 3. 

5.4.5 Gluing both normal and abnormal parts 

The tangent space to the sphere near the abnormal directions is described by 
the results of Section |5~T1 namely Theorem 15. 21 

Let A be the end-point of the abnormal trajectory and let K(r) be the first 
order Pontryagin's cone evaluated at A, r small enough. Let e i— ► a(e) be a C 1 
curve on the sphere S(0, r), a(0) = A, e ^ 0. Assume the following : 

1. a{e) c 5(0,r)\E r for s ^ 0. 

2. a(e) n L = 0, where L is the cut-locus for geodesies starting from 0. 

Then the tangent space to the sphere evaluated at a(e) tends to K(r) when 
e -> (see Fig. H5J). 




Figure 29: 



5.4.6 Lagrangian splitting and the Martinet sector 

Definition 5.5. We call Martinet sector of the Martinet sphere the trace of 
the ball -8(0, r) with the Martinet plane identified to y = 0. 

A precise description is obtained if we use the pendulum representation of 
Section 4, where the metric is truncated to order : g = (l + ay) 2 dx 2 + (l+(3x+ 
jy^dy 2 . The abnormal geodesic is strict if and only if a ^ 0. The pendulum 
equation is : 

6" + sin 6 + e/3 cos 8 9'+ e 2 a sm9(a cos 8 - j3 sin 8) = 

where e = ^= is a parameter. Cutting by y — induces a one-parameter 
section : 

S : 8' = e(acos0 + /3sin0) 
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The trace of the sphere with the Martinet plane near the end-point A = (— r, 0, 0) 
of the abnormal direction is described in Section 4. We take the first and second 
intersections of the pendulum trajectories with S, see Fig. I5U1 




Figure 30: (a ^ 0) 

In the conservative case, the curves D\, D 2 correspond to oscillating trajec- 
tories of the pendulum, and the curve C\ corresponds to rotating trajectories. 

Only one of the curves C\ , D 2 belongs to the sphere (this is D 2 on the figure) 
and their respective positions depend on the Gauss curvature of the restriction 
of the metric g to the plane (x, y). 

Contacts are the following. 

Proposition 5.14. Let Z — ^ and X = ^±r. Then : 

d,D 2 : Z = (± + 0(r))X 3 + o(X 3 ). 

D i : Z=-^X* + o{X*). 

This allows to describe the Martinet sector in the ball. 

Proposition 5.15. In the strict case a / the Martinet sector has the following 
properties. 

1. It is the image by the exponential mapping of a non compact subset of the 
cylinder : (9(0), A), A — * 00. 

2. It is homeomorphic to a conic sector centered on the abnormal line. 

3. It is foliated by leaves D\,E\ in the spheres 5(0, e), e ^ r, which glue 
according to Fig. 
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Figure 31: 




Figure 32: 



The Lagrangian splitting In the pendulum representation, the transport of 
T*^M n {H n = i} by the normal flow has the following basic interpretation, 
see Fig. E3 

The section splits into two parts S\ , S2 which represent the splitting of the 
fiber T*, Q sM into two Lagrangian manifolds. 

5.4.7 Microlocal invariants 

The pendulum has two singular points F — (0,0) and S — (0, 7r). The local 
analysis is as follows. 

• Near F, the linearized system is a focus whose eigenvalues are : 

e + • [71 2<p 2 7 

and is a perturbation of the linearized pendulum 8" + 9 = of the flat 
case. 

• Near S, the linearized system is a saddle whose eigenvalues are : 

V± = § ± ^/l + ^-a 2 ) = ±1 + y + o(e) 
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where e = -^=. It is a perturbation of the flat case r\± = ±1 which is 
resonant. 

In order to compute the sector we use the spectrum band r/± which is stable by 
perturbation. When we compute the sphere we have to compute an averaging. 
This is much more complex. 

The linear pendulum appears already in the contact case and the existence 
of the focus reflects the existence of a contact sector in the Martinet sphere. 

5.4.8 The sectors of the Martinet sphere 

In JHj was described the Martinet sphere in the integrable case by gluing sectors. 
We have three kinds of sectors : 

• A Riemannian sector R located near the equator, image of A = 0. 

• A contact sector around C, where C is a cut-point. 

• A Martinet sector around A, where A is an end-point of the abnormal 
line. 

The sectors are represented on Fig. 1331 




generic contact case generic Martinet case 

Figure 33: 

The microlocal invariants of the SR balls are the following : 

• Spectrum band corresponding to the focus F. 

• Spectrum band corresponding to the saddle S. 

• Invariants connected to the family of Riemannian structures a(q)dx 2 + 
b(q)dy 2 induced on the plane (x, y), where z is taken as a parameter. 

The conjugate points accumulate in the flat case along the abnormal direc- 
tion, see P]. 
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5.4.9 The n-dimensional case 

Our results except the precise asymptotics of Section 4 can be generalized to 
the n-dimensional case to define a Martinet sector in the SR ball. Indeed : 

• L 2 -compactness of SR minimizers (see £Q) allows to bound the number 
of oscillations of Lagrangian manifolds. It appears in our study by taking 
only the first and second return mapping to compute the sphere intersected 
with y = 0. 

• From , using a normal form, we can cut the SR ball by a 2-dimensional 
plane to identify a Martinet sector which splits into two curves : a curve 
Di obtained by using minimizing controls close to the reference abnormal 
control in L°°-topology ; a curve E± obtained by using controls close to 
the reference abnormal control in i 2 -topology, but not in L°° -topology 
(see Fig. 



The picture explains well the consequence of the existence of abnormal min- 
imizers in SR geometry. Contrarily to the classical case we cannot straight the 
geodesic flow near the abnormal direction to form a central field. 

6 Conclusion 

Our analysis explains the role of abnormal geodesies in SR geometry. It is based 
on the Martinet case. Using our gradated normal form of order 0, the geodesies 
foliation is projected onto a one-dimensional foliation in a plane which corre- 
sponds to a one-parameter family of pendulums. In this space the abnormal line 
projects on the singularities of the foliation. The computation of the sphere 
in the abnormal direction is related to the computation of return mappings 
evaluated along the separatrices of the pendulum. We have computed asymp- 
totics, using techniques similar to the ones used in the Hilbert's 16th problem. 




abnormal line 



Eh. 



Figure 34: 
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The computations are complex, because it is a singular perturbation analysis. 
In these computations one needs to consider geodesies C 1 -close to the abnor- 
mal reference trajectory on the one part, and geodesies which are C°-close, 
but not C 1 -close to the abnormal reference trajectory on the other part. Our 
asymptotics are not complete in the latter case and this requires further studies. 
Moreover the techniques have to be adapted to analyze the general case when 
the geodesies equations are not projectable. This leads to stability questions 
about our asymptotics. 

The projection of the geodesies flow onto a planar foliation, valid at or- 
der in the Martinet case, is useful to compute asymptotics but is not crucial 
from the geometric point of view, and Martinet geometry is representative of 
SR geometry with abnormal minimizcrs. The existence of such minimizers im- 
plies hyperbolicity seen in the pendulum representation as the behaviors of the 
geodesies near the separatrices. The general geometric framework to analyze 
SR geometry is Lagrangian manifolds. Here hyperbolicity due to the existence 
of abnormal directions is intcrpretated as a splitting of the Lagrangian fiber 
T* Q M when transported by the normal flow. To construct the sphere in the 
abnormal direction we must glue together the projections of several manifolds 
in the cotangent space. 

The link between the computations of asymptotics and Lagrangian geometry 
is the Jacobi fields which allow to compute contacts for the return mapping 
underlying our analysis. In general the evaluation of a return mapping in the 
analysis of a differential equation interpretated as a transport problem is original 
and source of further studies. 

The question of the category of the SR Martinet sphere is still open. The 
Martinet sector is homeomorphic to a locally convex cone and we have given a 
qualitative description of its singularities. In the integrable case the SR sphere 
is log-exp and this leads to a smooth stratification of the sphere. In general we 
conjecture that the sphere is not log-exp, belongs to some extended Il'Yashenko's 
category, and is still C 1 -stratifiablc. It is an important question connected 
to Hamilton- Jacobi equation which has many applications in physics (optics, 
quantum theory). For control theory, SR geometry is part of optimal control. 
Moreover our study is related to the stabilization problem. 
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